General Disclaimer 


One or more of the Following Statements may affect this Document 


• This document has been reproduced from the best copy furnished by the 
organizational source. It is being released in the interest of making available as 
much information as possible. 


• This document may contain data, which exceeds the sheet parameters. It was 
furnished in this condition by the organizational source and is the best copy 
available. 


• This document may contain tone-on-tone or color graphs, charts and/or pictures, 
which have been reproduced in black and white. 


• This document is paginated as submitted by the original source. 


• Portions of this document are not fully legible due to the historical nature of some 
of the material. However, it is the best reproduction available from the original 
submission. 


Produced by the NASA Center for Aerospace Information (CASI) 



FEB 20 1973 


NASA CR 120849 
BOEING DS8611-1 



ACOUSTIC RADIATION FROM LINED, 
UNFLANGED DUCTS - 

ACOUSTIC SOURCE DISTRIBUTION PROGRAM 

BY 

R. J. BECKEMEYER AND D. T. SAWOY 


DECEMBER 12, 1971 


the COMPANY 

WICHITA OiVtSION • WICHITA. KANSAS S72IO 


(N»Si-CB-12C8a9) ACOCSTIC RADIATIOS IHOH 
LIKED, UNFIANGED DUCTS: ACOUSTIC SCUECE 

DISTBIBUTICB PBOGEAE (Boeing Co., Hichita, 

Kans.) 202 p HC $7.25 CSCL 20A 

G3/71 


S75-28849 


Unclas 

29846 


PREPARED FOR 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 





NASA LEWIS RESEARCH CENTER 
CONTRACT NAS 3-14321 
H. BLOOMER. PROJECT MANAGER 



FOREWORD 

The work described herein wai done by The Boeing Company, 
Wichita Divifion, under NASA contract NAS 3-14321 with Mr. 
H. Bloomer, V/STOL and Nolle Division, NASA - Lewis 
Research Center, as Project Manager. 


ABSTRACT 


An acoustic radiation analysis has been developed to predict the far 
field characteristics of fan noise radiated from an acoustically lined 
unflanged duct. This analysis is comprised of three modular digital 
computer programs which together provide a capability of 
accounting for the impedance mismatch at the duct exit plane. This 
report discusses the Acoustic Source Distribution Program whose 
relationship with the other two modular reports of the analysis is 
illustrated on the following page. 

Admissible duct configurations include circular or annular, with or 
without an extended centerbody. This variation In duct 
configurations provides a capability of modeling inlet and fan duct 
noise radiation. 


DEVELOPMENT OF ACOUSTIC RADIATION ANALYSIS 
OF TURBOFAN NOISE FROM LINED, UNFLANGED DUCTS 

OVERALL REPORT ORGANIZATION 
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1.0 


INTRODUCTION 


1.1 Background 

This report discusses the Acoustic Source Distribution Program which forms the 
foundation for a package of modular computer programs which together provide a first 
generation capability for the analytical prediction of the far field noise being radiated 
from an unflanged inlet or exhaust duct. 

The development of these programs has its origin in the need for a method of predicting 
the far field characteristics of fan or compressor noise radiating from an acoustically 
treated duct. Among the techniques available in the literature are rigorous solutions for 
the case of sound radiation from unflanged, hard-walled circular pipes both with and 
without mean flow (uniform cylindrical jet) published by Levine and Schwinger, 
Reference 1; and Carrier, Reference 2 (valid only for plane wave propagation in the 
pipe); solutions for higher duct modes (no flow) by Vajnshtejn, Reference 3; and 
extensions to References 1, 2, and 3 by Lansing; Drischler and Pusey, Reference 4; 
where in the closed form solutions for higher duct modes (with flow) have been 
numerically evaluated. However, these solutions are not applicable to the more complex 
problems of noise radiated from lined and annular ducts. 

A method which can be applied to more general cases is the distributed source technique 
in which a surface distribution of elementary acoustic sources is determined, by 
numerical means, such that appropriate boundary conditions are satisfied on the surface. 
By distributing sources over a circular or an annular aperture in an infinite plane, 
Morfey, Referenc.e.s 5 and 6, and Tyler and Sofrin, Reference 7, were able to predict the 
acoustic radiation characteristics of flanged ducts. Researchers in underwater acoustics, 
for example, Chen and Schweikert, Reference 8, have modeled the acoustical radiation 
properties of a vibrating body of arbitrary shape by using distributions of acoustic 
sources over the entire outer surface of the body. Martenson and Liu, Reference 9, 
extended this work to the case of an engine nacelle. By distributing sources over the 
nacelle walls and the source (fan) plane, they represented both the propagating field in 
the nacelle and the radiating field exterior to the nacelle, thus, accounting for the 
diffraction effects at the lip of the nacelle exit plane. 

The Acoustic Source Distribution Program described herein is a further extension to 
these various source distribution methods which builds on the previous work to yield 
what is felt to be a basically more versatile method of approach. The program utilizes 
ducts of semi-infinite geometry, from which nacelles may be modeled by coupling of the 
inlet and exhaust radiation patterns in the far-field, Reference 10. The approach taken is 
to distribute acoustic sources over the duct walls and the duct exit plane. In this manner, 
edge diffraction effects are fully taken into account; however, the choice of the source 
and propagation (interior duct) model to be utilized is fairly arbitrary. The two problems 
are coupled at the exit plane through consideration of the duct termination impedance. 
Reference 11. 
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1.2 


Technical Approach 


The approach taken here involves the application of a numerical technique, the 
box-integration method, to determine a spatial distribution of acoustic sources over the 
duct walls and exit plane which satisfies the specified boundary conditions. The 
conditions imposed are, first, that the acoustic perturbation velocity normal to the rigid 
outer duct walls vanish, and second, that the velocity distribution over the face of the 
duct exit plane equals that associated with a given angular mode of sound propagation in 
the duct. 

Admissible duct configurations, shown in Figure 1, include unflanged circular or annular 
ducts, with or without an extended centerbody. Given the source distribution over the 
duct, the far-field directivity may be found by numerical integration of the Helmholtz 
integrals for the velocity potential at large distances from the duct exit plane, Reference 
10 . 


Note that this program assumes knowledge of the actual velocity distribution at the face 
of the duct. A given duct mode, when Incident on the duct exit plane, will generate a 
reflected wave composed of other duct modes because of the termination Impedance 
mismatch. Thus, the assumption that the velocity distribution at the exit plane is that 
due to an incident mode generated inside the duct is only a first approximation to the 
actual radiated noise. The Duct Termination Impedance Program, Reference 11, was 
written to provide a means for accounting for the impedance mismatch at the duct exit 
plane. 

There are several limitations to this first generation program. First, for the sake of 
expediency, mean flow within the duct has been neglected. It is anticipated that this 
simplification will not significantly affect results for noise propagation through a locally 
subsonic flow field of an inlet. However, radiation through an exhaust or Jet flow will 
undoubtedly not be prepresented properly by the present model. The second limitation 
to the program is restriction of the specified velocity distribution on the face of the duct 
to a Cos (m0) type of angular dependance. Note that coupling of different m modes in 
the far-field may still be accomplished by superposition of the source distributions, 
Reference 10. 

The assumption of a Cos {m^Jangular dependence does, however, restrict the Duct 
Termination Impedance Program to the extent that for a velocity distribution whose 
angular dependence Is specified by a combination of different m angular modes a 
separate analysis must be made for each set of angular modes. This limitation could be 
alleviated possibly by employing a collocation procedure in place of the box-method 
integration technique used in the present program. 

The most restricnve limitation to the program at present is one of size and corresponding 
computer time. As discussed in Section 4.2, the wavelength of the radiated sound affects 
both the maximum box size used in the program and the required integration distance 
along the duct wall. To adequately represent the spatial variation of the source strength 
distribution over the duct there should be at least two boxes per wavelength, {i.e., the 
box size decreases with increasing frequency). At the same time, however, in order that 
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the tar-field behavior at angles greater than 140 degrees be accuratelv rep it Is 

necessary that the integration on the duet outer wait extend back farther as the 
frequency increases. Together these two requirements indicate that at higher frequencies, 
a very large number of small boxes are required. 

A similar problem involves the requirernent for a very high number of angular boxes 
when attempting to represent the radiation characteristics of rotor-stator interaction 
modes at high angular mode numbers. 

As presently viewed, the primary change required in this first generation program is the 
need to develop a finite length duct represeriiAtion, This ehange is not a simple 
to the present program due to certajh inhereht rnathematicai problems^ Reference 12, 
but the basic logic and numerical techniques used presently are directly applicable to the 
anticipated formulation. 


MATHEMATICAL DEVELOPMENT 

Figure 2 illustrates the geometry of the problem. The outer duct walls and the duct exit 
face are each represented by distributions of mpnopoie singularities. Note that the 
source distributions along the duct walls are required to yield the unf|anged duct 
boundary condition of zero perturbation velocity normal to the outer duct wall. The 
velocity potential at an arbitrary point in the region outside the duct is given by: 

2 

Cr.d.rd) {!) 

where and represent the contributions of the source distributions over the 

duct exit plane, the outer duct wall, and the extended centerbody wall, respectively. 

Since ^ must satisfy the wave equation and represent outward propagating waves, it will 
have the form: 

0(R,t)= i 

R 


For discrete frequency noise, the velocity potential may then be written: 


Where 
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dS, = b d0,di, 
dSi ** a d^xdxj 


The ^j's are the unknown source strength distributions. 

0a — 00 

0 , = 

01 — 02 

The boundary conditions to be applied are: 

A. Zero normal velocity on the outer duct walls 
-d<f> 


U = 

u = 


dr 
dr 


Q 
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r=:b 

ziO 

r*a 

ziO 


(3) 


B. A prescribed velocity distribution over the face of the duct exit; 
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^ =V(r)GOs(m6!)e'‘^^ 


d r 
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Applying Equations (3) and (4) to Equation (2) yields the equations: 

dSj=V(r.S,0)e-'“' ® 
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METHOD OF SOLUTION 


The numerical technique chosen to solve Equations (5) through (7) is a 
"box-integration" method similar to that used in Reference 8. The method involves 
dividing the surface of the duct into a series of discrete areas and assuming that the 
source strength distribution over o given area or "box" is constant. The may then be 
taken out from under the integral sigh and the integrals evaluated. By reqUifing that 
Equations (5), (6) and (7) be satisfied at a finite number of control points, a set of 
simultaneous linear algebraic equations areqbtaihed for the unknown 

As the first step in solving the system of integral Equatiohs (5) through (7), the fact that 
the desired source strength distributions must be single valued |n angular dependence is 
utilized. Knowing the angular dependence of the modal velocity distributions in the duct 
tO beofthefOrm cos{m0},the ^iaheTaken as: ^ 

^ '■j’ (9) 

Thus by specifying the angular dependence of the the number of control points is 
vastly reduced since now Equations (6) through (7) need be evaluated at only one 
angular position. 


Next, divide the duct exit face into a series of M annular rings and the duct walls into N 
and N2 cylindrical segments as depicted in Figure 3. 


Considering for the moment only Equation (5), segmentation yields the relation: 

(r(,)cos{m0p) 5(r~ro 

•4 y»277_ ^ 
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where 


3(r-ro / 0-0o) 


Q if r#To or di^O^ 

1 if r= rj, and B=d^ 


Physical reasoning implies that 

lim 

Z-*^oo 


(111 


lim 0j(z) = O 
Z-fco . . 

since the contribution of the sources oh the duct exit plane to the vertical velocity on 
the duct wall at large distances becomes vanishingly small. This implies then that for 
practical purposes the integration boxes which extend to infinity may be neglected, that 
is, the integrations in the z directions may be truncated. Assume that the are constant 
over each of their respective boxes so that: 

M 

27r0g(I) cos(m0o)5 { r-rj,(9-f^o) 

T-1 

/ / f*^teerancl(z„^^,)Jcostn^^,')dS, (12) 

-f f [Integrand.(Zi,<9i!)Jcos(m^2!dS2— V(r) cos(m^)e ^ 

J«1 ’ *AZj *4=0 


For the sake of computational expediency, the control points for Equation (12) and the 
equivalent forms of Equations (6) and (7), are chosen at 0= 0, and at the center of the 
annular or cylindrical segments. Specifically, Equation (5) is satisfied atM points on the 
face of the duct. Equation (6) at N points on the outer duct wall, and Equation (7) at 
N2 points On the extended centerbody wall. By this procedure Equations (5), (6), and 
(7) become: 



where [D] is an M x M diagonal 

matrix whose elements are Djj = 2 ti; i 1, M, 

[Al is an IVl X (N + N2) matrix, 

[B] is an (N + N2) X M matrix, 

and [C] is an (N + N2) X (N + N2) matrix 

By expanding the partitioned matrix Equation (13), the unknown 0j's may be 
determined as: 
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( 14 ) 



A note on the method of integration used to obtain the elements of the coefficient 
matrices A, B, and C in Equation (14) is in order. 

First of all, the annular and cylindrical segments of the duct are further subdivided into 
L angular regions. The integrals are then evaluated over each segment by using an even 
order Gaussian quadrature scheme. The order of the quadrature is increased by 2 and the 
integral values compared. If agreement is not obtained, the order is again increased and 
the procedure repeated until the integral value converges. 

During development of the program difficulty in obtaining convergence was noted 
whenever the integration box was equal to or, in the case of small blocks, closely 
adjacent to a control box. This is a result of the singular nature of the integrals in 
Equations (5), (6) and (7) when Rj = 0, and is of concern since the number of arithmetic 
operations increases as the square of the order. For this reason, an approximate closed 
form solution for the integration on the outer duct walls for boxes centered on 0 = 
0*^ was derived. The procedure is valid only if these boxes have an angular width, 40 such 
that cos (m40) ~ 1, thus restricting the size of these boxes. Under this assumption, the 
closed form value contributed by a control block to its own vertical velocity on the duct 
outer wall (orcenterbody) is: 


400j for m - 0 

{ mAd \ 

2 sin \ 2 / (5 for m > 0 

m j 

Also, the contribution of a control block on the outer duct wall (centerbody) to the 
vertical velocity on another control block on the outer duct wall (centeibody) is found 
to vanish. Figure 4 shows the allowable variation of 40 with angular mode number for 
control boxes on the duct wail. 

The only other integrations which have difficulty in converging are those which 
determine the contribution of control blocks on the face of the duct, to the vertical 
velocity at control blocks on the duct walls (and vice versa). This problem of 
nonconvergence exists only for the control blocks on the face of the duct which are near 
the intersections of the exit plane and the outer and inner walls. As closed form 
evaluations for these integrals have not yet been obtained, the only solution here is to 
allow higher quadrature orders, and to accept the Increased computer time required for 
convergence. 
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4.0 


RESULTS 


4.1 Application 

The numerical techniques developed in Section 3 are applied to a discrete frequency 
radiation problem by the following procedure. First, a mathematical representation of 
the inlet or exhaust duct is determined by selecting one of the three admissible 
configurations shown in Figure 1. Next the duct segmentation parameters are selected, (a 
detailed discussion of the proper selection of box breakdown is given in Appendix 1). 
Finally, the velocity distribution on the duct exit plane must be specified (hard or soft 
wall) since the acoustic modes are specified by their velocity distribution. For example, 
to analyze a duct radiation problem of discrete frequency noise emanating from a lined 
duct, the soft wall modes and their respective contributions to the radiated sound energy 
must be determined. Far field directivity patterns are determined for each mode and 
then combined to yield the total noise signature. In the event both an inlet and fan duct 
are to be modeled, their respective far field directivity patterns can be calculated and 
then combined to determine the total noise signature of the nacelle. 

Several examples of applying the program to each of the three admissible duct 
configurations shown in Figure 1 are presented. The resulting source distributions and 
directivity patterns are used to evaluate the utility and accuracy of the program. For 
each example problem presented the frequency and outside duct radius were chosen to 
ensure that only the plane wave mode would propagate down the infinitely long duct 
Values of the acoustic and duct segmentation parameters are presented in Table 4.1. 


TABLE 4.1 

TABULATION OF ACOUSTIC AND DUCT SEGMENTATION PARAMETERS 


PROBLEM 

FREQUENCY 

PARAMETER 

kb 

RADIUS 

RATIO 

a/b 

SEGMENTATION PARAMETERS 
EXIT OUTER CENTERBODY 

PLANE WALL WALL 

INTEGRATION 

LIMIT 

1 

1.0 

0.0 

16 

40 

0 

4.1 

2 

1.0 

0.5 

16 

40 

0 

4.1 

3 

1.0 

O.S 

10 

10 

10 

1.375 

4 

1.0 

0.5 

10 

20 

20 

2,6 


Example Problem 1 

Noise radiating from a nacelle inlet Is represented by a circular duct configuration. The 
resulting source distributions on the duct exit plane0Q, and outer duct walI0^, are 
presented in Figures 5 and 9 respectively. To illustrate the behavior of ^ and ^-j, two 
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plots of each source distribution are presented. The smaller plot is a duplicate of the 
large one for an expanded scale. Note the extremely large slope of these distributions 
near the lip of the duct. This behavior which is due to the discontinuity in the velocity 
distribution, is discussed in detail in Appendix 1. Aside from the singularity near the 
duct lip, the behavior of the outer wall source distribution, (P-j, along the duct wall is 
important For a proper box breakdown, should exhibit a strong decay near the duct 
lip with small amplitude oscillations along the remainder of the duct. A normalized 
directivity pattern is shown on a polar plot in Figure 15. This pattern is compared with 
the results obtained by a numerical calculation based on the closed form solution given 
by Lansing, Reference 4. 

Example Problem 2 

The program was applied to an annular duct configuration to determine the effect of a 
centerbody on the far field directivity pattern. Source distributions for the duct exit 
plane and the outer duct wall are presented in Figures 6 and 10 respectively. The 
discontinuity in the source distribution on the duct exit plane at r/b = 0.5 is due to the 
discontinuity in the velocity at the inner duct wall. A comparison of the resulting 
normalized directivity pattern shown in Figure 16 and the normalized pattern of the 
circular duct shows the effect of the inner duct wall on the far field radiation. 

Example Problem 3 

Noise radiating from a nacelle fan duct Is modeled by extending the inner centerbody 
past the duct exit plane. Source distributions for the duct exit plane, outer duct wall and 
extended centerbody wall are shown in Figures 7, 11, and 13 respectively. Singularities 
in the exit plane source distribution occur near the extended centerbody wall as well as 
near the outer lip of the duct. The extended centerbody source distribution ^2 exhibits 
several large oscillations before decaying. The effect of the extended centerbody wall on 
the normalized directivity pattern is shown in Figure 17 to modify the circular duct 
pattern near the 90° region. 

Example Problem 4 

The annular duct with an extended centerbody is modeled with a different box 
breakdown than that used in Problem 3. The box breakdown of the extended 
centerbody wall was refined near the duct exit plane and additional boxes on the 
centerbody carried the integration farther back on the wall. Source distributions for the 
duct exit plane, outer duct wall and extended centerbody wall are shown in Figures 8, 
12 and 14. The directivity pattern of Figure 18 shows the effect of modifying the box 
breakdown on the duct. Finer breakdown yields a pattern that is closer to circular duct 
pattern than was observed in the coarse breakdown of Example Problem 3. 


4.2 Limitations 

There are several limitations to this program. First, mean flow within the duct has been 
neglected. Although application of this analysis to a nacelle radiation problem does not 
yield an exact model, it is anticipated that this simplification will not significantly affect 
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results for noise propagation through a locally subsonic flow field of an inlet. However, 
radiation through an exhaust or Jet flow will not be represented properly by ttie present 
analysis. 

The second limitation to the program is restriction of the specified velocity distribution 
on the face of the duct to a Cos (m9) type of angular dependence. This assumption does 
not affect the capability of the program to combine different m angular modes in the far 
field by superposition of the source distributions, Reference 10. It does, however, 
restrict the Duct Termination Impedance Program, Reference 1 1 to the extent that for a 
velocity distribution whose angular dependence is specified by a combination of 
different m angular modes a separate analysis must be made for each set of angular 
modes. 

The most restrictive limitation to the program at present is one of required storage size 
and corresponding computer time. Both of these digital computer usage parameters are 
related to the wavelength of the radiated sound waves. To adequately represent the 
spatial variation of the source strength distribution over the duct, there should be at least 
two boxes per wavelengtii (6 boxes per wavelength is more desirable). In order that the 
far-field behavior at angles greater than 140 degrees be adequately represented, it is 
necessary that the integration on the duct outer wall extend back farther as frequency 
increases. Together these two requirements indicate that at higher frequencies, a very 
large number of small boxes are required. 

Time and storage requirements depend on the number of integrations required to fill the 
influence matrix in Equation (13). The number of integrations varies with the square of 
the number of boxes and computer run time is dependent on the number of integrations. 
When convergence of the numerical quadratures becomes a problem the program 
requires a considerably larger amount of time to compute the matrix elements. Core 
storage required by the computer is equal to twice the square of total number of boxes, 
since the matrix elements are complex numbers. The number of boxes required on the 
exit plane of the duct becomes critical when radiation patterns of high radial modes are 
desired. At least two boxes per zero crossing of the velocity distribution is required, {i.e. 
the 6th radial mode requires at least 12 boxes. A similar problem involves the 
requirement for a very high number of angular boxes when applying the program to the 
problem of high angular modes. This requirement increases the computer run time since 
the number of numerical quadratures is greater. 
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5.0 


CONCLUSIONS AND RECOMMENDATIONS 


An acoustic radiation analysis has been developed that predicts the angular distribution 
of acoustic energy in the far field that is radiated from an unflanged duct. This analysis 
represents a rather unique analytical tool in that it can be applied to hard wall ducts 
(unlined) with an annular cross-section as well as circular. Both circular and annular duct 
configurations with their interior wails lined with noise suppression material can be 
analyzed. For all duct configurations with either hard or soft walls, the far field radiation 
patterns are corrected for the impedance mismatch at the duct exit plane. Also, of 
primary importance is the fact that the present analysis permits combination of inlet and 
exhaust radiation patterns in the far-field. 

In spite of its many advantages, this acoustic radiation analysis is limited in its range of 
application. For this reason, this analysis has been relegated to the role of a first 
generation analysis in a commitment to develop an analytical package which will predict 
radiation characteristics of a turbofan jet engine. The present version of the analysis 
reaches optimum efficiency when it is applied to low frequency noise radiating from an 
infinite duct. This and other limitations due to a Cos (m0) type of angular dependence of 
the velocity distribution and the required storage size and corresponding computer time 
are an inherent property of the "box method" numerical technique. 

A second generation version of the program could avoid these limitations by utilizing a 
collocation technique of evaluating the integrals. The primary feature of this technique is 
that a selected analytical expression is used to represent the source distribution. The 
time saving feature of the collocation technique is that it requires fewer evaluations of 
the integrals. The collocation procedure would enable the analysis to handle any 
combination of angular modes, a result which would remove the angular modal coupling 
limitation for the Impedance Program, Reference 11, and eliminate the need fora large 
number of angular boxes. 

In addition to the collocation procedure, a second generation version of the analysis 
should contain a method of accounting for a mean flow withm the duct. This method 
should account for the effects of velocity gradients on the distortion of the radiation 
field by the jet. 

A further refinement in the overall development of an analysis to predict turbofan noise 
would be the incorporation of a finite length duct representation. This change is not a 
simple extension to the present program due to certain inherent mathematical problems 
that involve the nonexistence of solutions of the integral equation formulation. As 
discussed In Reference 12, these problems can be avoided by using an alternate 
formulation. 

A finite duct length representation would require additional analytical considerations. 
The presence of higher radial modes generated at the duct termination and reflected 
back toward the source couid not be ignored. Their reaction with the noise source and 
subsequent propagation down the duct towards the exit plane requires that the problem 
must be analyzed by a systems analysis approach. 


‘MX ' 
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NOMENCLATURE 


Computer Program Variables 


A 

AA(U) 

B 

C(l) 

C2(l) 

CHK 

DD(I,J) 

D(l) 

INK 

IV 

K 

L 

MW 

M 

N 

N2 

NNZ 

PHKI) 

R 

RO,RI 

S 

TT 

Tl 

VII) 

W 

WDS 

XANG 

zz 


Centerbody Radius 
Matrix 

Outer Duct Radius 
Matrix 

Dimensionless Size of Ith Box on Duct Outer Wall 
Dimensionless Size of ith Box on Extended Centerbody Wall 
Convergence Criterion 
Matrix 

Dimensionless Size of Ith Box on Duct Exit Face 

Configuration Parameter 

No. of Velocity Distributions to be Read In 

No. of Polar Angles at which Far Field Data is Required 

No. of Angular Boxes 

Angular Mode No. 

No. of Annular Boxes 

No. of Cylindrical Boxes on Duct Wall 

No. of Cylindrical Boxes on Extended Centerbody Wall 

No. of Velocities to be Computed by VE LI NT 

Source Strength Distribution 

Radial Coordinate 

Distance Between Point on Duct and Field Point 
Sound Speed in Air 
Angular Coordinate 

Angular Coordinate Integration Variable 
Velocity Distribution on Exit Face 
Frequency Radians/Second 
Wavenumber W/S 

Nondimensional Angular Size for Control Box 
Axial Coordinate 
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INPUT DATA 


Since the accuracy of the source strength distribution obtained from this program is critically 
dependent on the method of segmentation of the duct walls and exit face, a brief discussion of the 
pertinent points is in order. In addition to the requirements that the control boxes have a small 
angular width such that Cos(m40)»1.O and that there be two or more boxes per wavelength of the 
radiated sound, there are restrictions on the size of the annular and cylindrical su>gment heights (Ar 
and Az, respectively) in the region of the intersection of the duct walls and the exit plane. This 
requirement must be satisfied, primarily, to accurately compute the near field pressure (over the 
duct exit face). 

The discontinuous velocity along this intersection results in the 0:'s being singular there. (An inner 
asymptotic expansion of the wave equation near the intersection has revealed that the singularity is 
one in slope of the 0j versus z or r curves. Specifically, as the intersection is approached, the 0j 
behaveasjx|ln|x|, x going to zero, where in denotes the natural logarithm.) 

Adequate representation of this singular behaviour requires a fine breakdown in the region of the 
intersection. For a uuct of radius b, the breakdown near the intersection should consist of several 
boxes ranging in size from 0.01b up to 0.05b. For all duct configurations the fine breakdown /vould 
be necessary near the duct lip on both the face and the wall. For an annular duct it is also necassary 
on the face on both sides of the inner annulus wall, whereas for the annular duct with an ex<;ended 
centerbody, the fine segmentation is necessary on the face and centerbody wall near their inte "section. 

As a final note, the box sizes on the face in both the radial and angular directions must naturally 
decrease in size when higher order modes with many modal lines are represented. 

CARD OR 

CARD SET FORMAT 

1 18A4 

2 NAMELIST 


3 3012 


(0 - Circular or Annular Duct, 

1 - Annular Duct with Extended Centerbody) 


DATA 

Title of run 

Logical Parameters: PUNCH = T(F) yields (does not yield) punched 
card output to be used in the directivity index program. 

OLD = T(F) bypasses (does not bypass) the integration loops and 
uses (does not use) previously generated matrices in the MATX 
program to compute PHI distributions for new velocities. Delete card 
(sets) 3 through 7 if OLD = T. 

Integer Parameters: 

L = No. of Angular Segements 
M = No. of Annular Segments 
N = No. of Cylindrical Segments 
NN = Order of Initial Gauss Quadrature (NN Even) 

LP = No. of integrations in Integration Routine 
MW = Angular Mode Number 
INK = Configuration Parameter 
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4 7F10.3 


5 7F10.3 
(Card Set- 

N Values, 

7 Per Card) 

6 7F10.3 
(Care! Sel- 

T42 Values, 

7 Per Card) 

7 7F10.3 
(Card Set- 

M Values, 

7 Per Card) 


8 212 


9 6E12.5 

(Card Set- 
M Complex No's. 

3 Per Card 


K = No. of Angles at which Far Field Pressure Values are to 
be computed. 

N2 = No. of Cylindrical Segments on Centerbody. 

A = Centerbody Radius 

(A = 0.0 for Circular ^ Annular Ducts) 

B = Outer Radius of Main Duct 
W = Frequency, Radians/Second 
S = Sound Speed in Units Consistent with those of A & B 
CHK = Convergence Criterion for Numerical Integration 
XANG - Nondimensional Angular Size for Control Box 
A6 Radians = 27 t (XANG) (See Figure 4). 

C(l), I = 1, N, the size of the cylindrical segments on the Main Duct 
Wall in the Z direction, Nendimensionalized with respect to B, 
Numbered with C(1) closest to the duct exit plane and I increasing in 
the negative Z direction. 

C2(l), I = 1, N2 In put this set only if INK = 1. Size of cylindrical 
segments on centerbody wall in Z direction, nondimensionalized 
with respect to B, numbered with C2(1) closest to the duct exit 
plane and I increasing in the positive Z direction. 

D(l), I = 1 , M, the size of the annular segments on the duct exit plane 
in the radial direction, nondimensionalized with respect to B, 
Numbered D (1) at the duct centerline (or adjacent to the 
centerbody wall-duct exit plane intersection) and increasing to D(M) 
at the lip of the duct. 

IV 

IV 0, IV = No. of velocity distributions to be read from cards. IV = 
0, velocity distributions to be computed as hard-wall modal 
distributions from VELINT. 

V(l), I = 1, M, Velocity* distribution on face of duct 

V(K) = velocity corresponding to segment D(K). - In put this card set 

only if rv>0. 


10 3612 NNZ - No. of Velocity Distribution to be computed by VELINT. 

In put Card 10 only if IV = 0. 


•NOTE: For annular duct without extended centerbody, the velocity for the boxes on the inner 
annular area must be set equal to zero. 
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PROGRAM FLOW 


Flow charts for the box method program and associated subroutines are given in Appendix 2. 

The main program consists of the calling program and subroutines CAL200, CAL300 and CAL7Q0. 
Here the input data are read in and the box method integrations performed. Upon completion of 
the integration, subroutine MATX is called. This subroutine carries out the matrix manipulations 
required to solve for the source strength distributions. 

Subroutine GAUSS Is used In the integration and consists of tables of Gaussian weights and 
abscissas. Subroutine INVT performs the inversion of complex matrices. Subroutine VELIISIT is an 
optional routine which computes modal velocity distributions for hard walled annular ducts. Thus, 
the velocity distributions may be input directly or computed by VELINT. 


OUTPUT 

Output quantities consist of the following: 

A. input data printout 

B. Results of the integration program, consisting of the contribution of each box to the 
vertical velocity at each control point. 

C. Results of the matrix operations during which the 0j's are computed. (Optional) 

D. The values of the source strengths at each box. 


EXAMPLE PROBLEMS 

Appendix 8 illustrates the input preparation and output format of the program. A listing of the 
input data are given for each example problem discussed in Section 4.1. Complete printout (with 
the exception of the operations performed by MATX) is shown for example Problem 1. 
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PROGRAM FLOW CHARTS 
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15 


lt-/0) 

I 

( 

1 

I 

I 

I 




16 


I Rtn p I 

l RARCtnftBB/Z.DO ft I 
I AINK I 

« — . — - , ...■■■I - m 

1 

I 


L70 


1 BESHtn 


23 

O.DO 1 


I 

H- 


50 


ft NO 

END OF DO ft-ft 
LQOP ^ ft I 

1 


* 

... 

IY£S 

. 45 

1 

. 15 

I 

I 

1 

*« . 

1 

I 

1 

Z5 


f HRITE^TO DEV / 
/ VIA FORMAT 
/FROM the list / 


I NOTE 26 

ft LIST X H« ft 

I 

1 NOTE 27 

««*ftft*ft«»*ft 

ft OEGIN DO LOOP ft 
ft 54 J * I, KARC ^ ft 


52 




17 


I IE p 1 - I I 

1 Rt n a Rl IE) ft I 
f (RAROtlEJ ft 1 

I RARG{in«BS/2.D0 1 

» ■ ' — ■ — X 

I 

1 

- - tils CANNOT 

BE USED V11H CIRCULAR 
DUCT CATA 


1 

■M 


IB 

ft 


*1NK ,EQ. 0 ft TRUE 
.ANO, rFiI 


nnu* X J 

•LT« AA 


54 


-i! 

1 


/ WRITE TO OEV < 
/ 6 / 
t VIA FORMAT f 
t 53 / 

i FROM THE_IIST / 

"l 
I 


LIST 

RIJ) 


42.03- 

II 


f NOTE 2 

a AARGIJ) » 
j beshjJ ^ 

J 

J 

1 


END OF DO 
LOOP 


IVES 
I 
t 
I 
t 

-n 


ifalse 

I 

1 

I 


NOTE 3L 


CONTINUE 


BARC ft 

RC n«XRMN/liD 


19 


t 

I 

I 32 
ft EXIT 


•':SND!"rN8** *i:hl 

ft .Eg. 01 • 


-‘T 


BESJ ft l.DO 
BEST « O.DO 


33 

I 

I 

i 


FALSE 


14 { 
tS I 

I. I lEARCtnWTBESJt H 
10 I BESYtBiiH.eunni h 


^CBESJY 


rtflUHtB 

IEHR> 


1 NOTE 34 


CONTINUE 


I BESMin a BESJ 
I qHH»8ES V 

1 

I 


35 


82 


/ OQESJY / 


/ 


2DH915 as - 29 


— — a I 

} -i.- "i 

I l^l9926939e97930C | 

I 
I 
I 


* « TRUE 

X *GT. O.DO •- 



. - i 

04 

* 

EOKTINOE 

« 

• * « 

9*9 999 

m 9 


I 



I 

05 

I 

tERR « 1 

1 

I 


I 

^ 

FAC - US 

1 


FALSE • 


X «GT# 99.0 


IFALSE 

1 

t 

t 

I 

I 

1 


03 
• * 
I 


.4B.26. 

• So 


ITFOE 

J 
1 
I 
1 



( 

07 

! 

FAC - -5 

i 


-u‘ 


* 


08 


I NX * 3.9 f 


Jt2*N>C ♦ 
0»00 


FAC*X 
2i « 


J(2*NX 

UO - 


' II 
29 


1 

t 

1 

c « o*co 

JO X O.DO 

1 

\ 


I 

I 

1 

09 


Jl « 0.00 

I 


YO - O.DO 

I 


n - O.DO 

I 






, HDTE 

10 

9 


9 

9 

BECZN DO LOOP 

9 

9 

9 

ZO NB 9 L« NX 

-li 

9 



1 L 


N - NX » 1 - NB 

1 

1 

RN « OFLOATINI 

1 


J124N» » 

t 


2.00»[Z.DC«RN 9 

I 


UQ0]»J'2BN f 

1 


I) /X • JU9H 9 2) 

I 


5 

12 


Ji2*H ~ \i » 
4.00* RN*J I2*N} 

I 

I 


/% *- J12*N ♦ IJ 

I 



I 


YO ^ YC , 
l-U»*N*J<2»hJ/RN 



1 



I 

I 

13 

1 

n « ri ♦ 

I 

1 

t-n 

1 

I 

«*N*I2.Q0»RN 

* t 

I 

. . UDO>»J< 2*N 

♦ t 

1 

tl /Ith/U.DO * 

RN| I 


I 

20 1 

u 

i 

K- 

c * c ♦ 

Z.OO^JI 2«H] 

1 

1 


END DF 00 
^ LOOP 


res 


[ i 


OEIGINAn PAGS'li 
OF POOB QUAiM, 


JDO • 

2«00*JCli/X - 
JC2I 

C • JQO » C 
YO • 

(PLGGlX/2.00] * 

S772I566A9O153D0I 
»JOD - 2. P0»YD 

X 

I 

1 


lb 

— 9 

t 

1 

I 

1 

I 

I 

1 

1 

I 

I 

I 


Y1 ■ 

tOLOGU/2*QOJ - 

9227fi43350<’,8f700) 
•Jl IJ - JOD/X - 
Y1 


VO 

VI 


I 
t 
1 
1 

2.DC«Y0/C/PI I 

2,U0yi/C/Pl T 

1 
t 
Z 


JO « JQO/C 

vii» ■ ri 


1 VtZJ - 1 

I 2.D0«rt LWX - VO I 


NhX X 2 » NO 


1 NOTE 19 

«**«»*9*» 

BECIN 00 LOOP * 
30 I X l, HKX 



1 

» 9 


^ I 

20 

T 

YU 2) - 

1 

I 

2.DQ«0FL0ATt 1 ♦ 

1 

1 

IJ«TCI *.l»/X - 

t 

k 

Y[ II 

1 


I 

30 ! 

21 

I 

Jt n * jui/c 

1 

»■ 


Z2 


IVES 

] 


23 


9 NO .EO. O 


IFALSE 


I 


BJ - JINO) 
BY ■ VtNO> 
NP ■ NO « 1 
BJJ > JtNPl 
BYT ■ VINPI 


2^ 

.tt 

1 

[ 

I 

I 

t 

I 


APPENDIX 3 

MAIN PROGRAM LISTING 


PRECEDING PAGE BLANK NOT HLMEO 
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C BOX METHOD IN'IEGRATION FOR SOURCE DISTRIBUTION OVER DUCT 

C ZDM9L5 f.1 - 61 

LOGICAL * 1 PUNCH/T/ ,OLD/F/ 

INTEGER *4 TITLC/*-Z»/r TITLEO/*R'/ ,DRECTP/02/ ,TITLE(18) 
1 ,SAVTP/03/ ,TITLC2/*+ZV ,MAXD/60/ 

REAL *4 ZERO/0.0/ 

DIMENSION C140J, 0(401 ,BR(60t60) 


,BI (60,60 I 


ASQ ,DTl, 0T2 


1 ,E(40) ,C2t40) 

COMMON /GAUS/ 

1 0SG1, WDS, ZL, RL, 

2 , ICN 

COMMON /MATGUS/ B2t AINK, PI 2, 

DUMMY IN SO COMMON'S WILL AGREE SIZE WISE 
N, H, MW, LP, L, NIM,N2 
BR,ai 


C, D, A, Br CHK ,C2 


DIRECTION NORMALIZED*, 




- * — * 


I flNK, 

COMMON /SAVl/ 

100 F0RMAT(30I2) 

101 FORMAT! 7F10. 31 
103 FORMAT ( 18A41 

106 FORMAT (*0*, 5X, 'BLOCK SIZE IN THE *,A?, 

1 ' WITH RESPECT TO B*' /(3X,4El5,4n 

150 F0RMAT11H1,5X,«B0X METHOD INTEGRATION FOR SOURCE DISTRIBUTION OVER 

1 DUCT',// 6X, 18A4 

2 // 6X, 'DUCT OR = ',F10.4,* FT •,3X,*DUCT IR =',F10.4,« FT 
3*,/,5X, 'FREQUENCY =',F10.4,* RAD/SEC 5 X ,' SOUND SPEED =',F10.4,2X, 
4'FT/S€C',/,5X, 'MODE NUMBER =',I3,//1 

NAMELIST /PRNTl/ L,M,N,NN,LP 
NAMELIST/ PRNT2/ N2, INK, CHK, XANG 

NAMELIST /RFAD/ PUNCH , OLD 

READ (5,103 ,END=2) TITLE 
OLD = .FALSE. 

READ ( 5, READ) 

IF (OLD} 

I READ(5,100 ,END=2) L, 

IF (L ,LE. 0) 

IF (M-t-N4-N2 pGT. MAXD) 

READ(5,lOl) A,0,H,S,CHK 
WRITE(6,l5ni TITLE, B 
WRITE( 6, PRNTl) 

WRITE! 6, PRNT2) 

READ(5,10U IC(I), 1 = 1, N) 

WRITE (6,106) TITLC , ( C ( I ) , I = 1 , N 1 
IF (INK .LT. 1) 

R FAD( 5,101 1 IC2 ( I ) , I =l,N2) 

WRITE (6,106) 

3197 READ! 5, 1011 (0(1), 1=1, Ml 

WRITE (6,106) TITLED , ( D( I 1 ,1 =1 ,H) 


M, N, NN, LP, MW, 


GO TO 75 
INK, K ,N2 

GO TO 2 

GO TO rooo 


,XANG 
A , W , S , 


MW 


GO TO 3197 


TITLC2, (C2{I 1 ,I=1,N2) 


INITIAL CALCULATIONS 
PT = 3. 14159 26 53 58979 
E( 11 = XANG 


/ a-1) 


* PI / tL-n 


f3T2 = (l.OEO - XANG) 

DO 50 I=2»L 
50 E( 1} = 0T2 

DTI = XANG ♦ PI « 2,0 
DT2 = 2.0^ (l.OEO - XANG} 

ICN=0 

PI2 = PI 2 
B2 = fi / 2.0 
BSQ = B * B 
WDS = W / S 
ASQ = A * A. 

AINK = A * INK 
CALL CAL 200 
CALL CAL 300 

IF (INK ,GF. 1) CALL CAL700 

; - - PAD WITH ZEROES SO THEY MAY BE READ IN REAL *8 

7000 WRITE (SAVTP) INK, N, H, HW, LP, L, NN, K , N2 

1 ,B2,ZERfJ, AINK, ZERO, PI2,ZERQ, A, ZERO, B,ZERO, CHK,2ER0 

2 ,W,ZFRQ, SfZERO, ( C( H ,ZERO , I =L ,NI 

3 ,{Dm,ZFRD, 1 = 1 , M) , (Ein , ZERO, 1 = 1, L) , t C2 ( II , ZERO, 1= I , N2 ) 
NUM = N H + N2 

WRITE (SAVTP) ((BR(I,J ),ZERO, BI ( I , J J , ZERO, I = l , NUM) , J= I , NUM » 

75 CALL MATX ( ORFCTP , PUNCH , SAVTP, TITLE, UNI TEE) 

GO TO I 
2 CONTINUE 
STOP 
END 
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SUBROUTINE CAL 200 


C ZDM915 62 - 61 

DIMENSION C(40),D(40) ,BR 1 60* 60 ) ,BTR ( 60) ,BI( 60 , 60 ) » BT I ( 60 I 


1 

fX!16), 

MT!16) ,C2!40) 


COMMON 

/GAUS/ 



1 

hSO 

, WDS, ZL, RL, ASQ ,DT1, 

DT2 

2 


, ICN 


COMMON 

/MATGUS/ 

B2, AINK, PI2, 

C, D, A« B, CHK ,C2 

I , 

INK* N, M, 

MW, LP, L, NN *N2 


COMMON 

/SAVl/ 

BR,B I 


- * 

- ♦ - 



FMW = 

FLOAT!MW ) 



DO 200 

JC=l,M 



IF! JC- 

1)210,210, 

2U 



210 R=DUCI* B2 + AINK 
GO TO 212 

211 JE=JC-1 

R=R+(0< JE) +DUCn*B2 

212 WRITE! 6, 153) JC,R 

153 FORMAT! 1H1»5X, ’BLOCK NO* , I 3, 5X, ’RADI US = •*F10-^,//) 
ICN=ICN+1 

IDN = 0 

CALCULATION! OF THE D MATRIX ELEMENTS 

on 220 JJ=I,M 
IDN=IDN+1 

IF! JJ-JC J221*?22,221 

222 BR( ICN* ION)= P 12 
GO TO 223 

221 BR! ICN, ION )=0.0 

223 BI! ICN, IDN 1-0.0 

IFIMOO! IDN, 45) .FQ.O) WRI TE ( 6 , 3718) 

3718 FORMAT! IHl ) 

WRITE! 6, 15 4) IDN,BR( ICN, ION) ,BI|ICN,I0N) 

154 FORMAHIH , 5X, 13 , 5X* F 1 0. 4, 2X , • + * ,2 X,F1 0. 4 ,2 X, * I * ) 

220 CONTINUE 

CALCULATION OF THF A1 MATRIX ELEMENTS 

DO 230 K = 1,N 
R2C = B2 * C(< ) 

IPIK-D231, 231,232 

231 ZZ=-B2C 
GO TO 23 3 

232 KL=K-1 

ZZ=7Z-!C£KL)+C(K))4B2 

233 IDN=IDN+1 

ZR = ABS !ZZ) 

BR! ICN, IDN )=0.0 
BI! TCN, IDN )=0.0 
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Z3^t 235, 237 


DO 270 1 = 1, L 
IF ( I - 2) 


234 

DT 


DTI 





TT 


0.0 

GO 

TO 

236 

235 

DT 

= 

DT2 





TT 

s 

(DTI + DT ) * .5 

GO 

TO 

236 

237 

DT 

= 

DT2 





TT 

= 

TT + DT 





236 B0TD4 = B * DT / 4.0 
DTD2 = DT / 2.0 
BDTD4C = RDTD4 ♦ C(K ) 

DO 250 J=l,LP 
ILK=J 

ILL=NN + 2*( J-U 
BTRIJ)=0.0 
BTII J)=0-0 

CALL GAUSSlXrWT, ILL) 

SANR^O.O 

SAM 1 = 0 .0 

DO 255 I r=l, ILL 

ll=ZZ- X (I n ♦ B2C 

7ISQ = ZI ♦ 21 

DABSZI = ABS (ZI) 

DO 255 JJ=1, ILL 
TI=rr- X{JJ) * DTD2 

RI = [(R-B * C0smn*^2 + BSQ ♦ ( SIN(n>)*+2 *■ ZrSQ )**n.5 

ART=R I^WDS 

CART = cast ART) 

SART = SINIA-T) 

SINT=Z I*CQSII-'^W«Tn / (RI*RI) 

DUMY = B * SIMT * HTIII) WT(JJ) 

SANR= -DUMY* ICART/RU SART*HDS) +SANR 

SANI=OUMY*(WOS*CART -SART/RI) +SANI 

255 CONTIMUE 

BTR( J)=SANR*B0TD4C 
BTl( J 1=SAMI*BDTD4C 

IF { J .!. T. 2 ) GO TO 2 50 

JL=J-l 

ABT = ABS( AHSIBTRtJ}) - ABSIBTRUU)) 

IF (ABT ,GT. CHK) GO TO 257 

ABC = ABSI ABS(BTUJJ) - ABStBTKJL))) 

IF (ABC .LE. CHK) GO TO 260 

257 IF( J-LP)250,253, 260 

253 WRITE (6,6500) ABT, ABC, TT, ZZ ,BTR(J), BTKJ) 

6500 FQRMAT(* IMTFGRATIOM DID NOT CONVERGE - RFAL DIFFERENCE =*, 

1 E13.5,/* IMAGINARY DIFFERENCE =', E13.5 ,5X* THET A =*,E13.5/ 

2 /» Z =• ,F13. 5, SX'BTR = • ,E13.5 ,5X* BTI =*,E13-5) 

250 CONTINUF 

260 KIL=ILK-l 


o n o 


BR( ICNt IDN)=BR( ICN,ION)+BTR<ILK) 

BI( ICN, { ICNi LON j-i-BTinUK) 

270 CONTINUE 

IF(MOO( ION, 45) .EQ.O) HR I TE (6,37181 
WRITE( 6, 15 5) IDN ,BR{ ICN, IDN) ,BI ( ICN, ION) , ZZ 
155 FQRMATtlH , 5X, 1 3 , 5X, FIO. 4, 2X , » + • ,2 X ,F 10. 4 , 2X , * 1 • ,8 X , » Z = •,F10.4) 
230 CONTINUE 

CALCULATION OF THE A2 MATRIX ELEMENTS 

IFI INK-U 200,500, 500 
500 DO 530 K=l,N2 
B2C = B2 *C2(K ) 

TFJK-U531, 531,532 

531 ZZ=B2C 

GO TO 533 

532 KL=K-l 

ZZ=ZZ<'IC21KL) + C2(KM * B2 

533 IDN=IRN+1 

ZB = ABS (ZZ) 

BR( ICN, ION )=0.0 
BI(ICN,IDN)=0.0 
DO 570 1=1, L 



IF 

( 1 - 

2) 

534, 535, 

534 

DT = 

DTI 




TT = 

0.0 


GO TO 536 

535 

DT = 

DT2 




TT = 

( DT] 

+ DT) * ,5 

GO TO 536 

537 

DT = 

DT2 




TT = 

TT 4 - 

DT 


536 

BDTD4 = B 

♦ DT / 4.0 



UTD2 

= DT 

/ 2.0 



BDTD4C = BDTD4 *C2(K) 

DO 550 J=1,LP 
ILK=J 

ILL=NN4-2*( J-U 
BTR(d )=0.0 
BT I( j )=n.o 

CALL GAUSS! X,HTi ILL) 

SANR=0,0 

SANI=0.0 

DO 555 I I- 1, ILL 

ZI=ZZ- Xtll) ♦ B2C 

ZISQ = ZI * ZI 

OABSZI = ABS (ZD 

ZIA = It * \ 

DO 555 JJ=l,ILL 
TI=TT- X(JJ) *0TD2 

RI = i(R - A* CDS(TI))**2 *■ ASQ ♦ I SIN(TI))*^2 ♦ ZISQ )**0.S 


90 


555 


557 

553 

550 

560 


570 


5 30 
?00 


ART=RI*W0S 
CART= COS(ART) 

SART= SIN(ART> 

SINT=ZIA * CnS{FMW*Tn / iRI^RI) 

DUMY = SIMT ♦ WTim * WTfJJl 

SANR= -DUMY«= (CART/RH- SART^WDS) 

SANI= (WDS* CART -SART/RIJ *OUMY i-SANI 

CONTINUE 

BTRlJ)=SANR=«‘a0T04C 
BTIl J)=SANI=«‘BDTD4C 

IF (J ,LT. 2 ) GO TO 553 

JL=J-1 

ABT = ABS( ABSI BTR( J M - 
IF ( ABT ,GT. CHK ) 

ABC = ABSI ABSI BTK J I) - 
IF (ABC .LE. CHK) 

IF(J-LP)550»553, 560 
WRITF 16,6500) ABT, ABC, TT, ZZ ,BTPU), BTI(J) 
CONTINUE 
KIL=ILK-1 

BR t ICN, IDN ) = BR( ICN, I DN l + BTR IILK) 

DU ICN, IDN )=BI (ICN, IDN)+RTI ( ILK) 

CONTINUE 

IF (MODI I ON, 45) .FO.O) WR I TE ( 6 , 3718 1 

WRITFl 6, 15 5) lONtBRl ICN, IDN) , R 1 1 ICN, I DN) , Z 7 

CONTINUE 

CONTINUE 

RETURN 

FMD 


ABSIBTRI JL) ) ) 

GO TO 557 

ABSIBTl IJU I ) 

G ) 560 


+ SANR 
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SUBROUTINE CAl300 


C ZDM915 63 - 61 

DIMENSION C(40),D(40I , BR ( 60 » 601 »BTR { 60) , BI ( 60 , 60 ) , BT I ( 60 J 


I 

,XC 16) , 

WTI16) ,C2!40) 



COMMON 

/GAUS/ 




1 

BSQ 

, WD S , ZL, RL , 

ASQ ,DT1, 

DT2 

2 


, ICN 



COMMON 

/HATGUS/ 

62, 

AINK, PI2, 

C, D, A, B, CHK ,C2 

1 » 

INK, N, M, 

MW, LP, L, NN 

,N2 


COMMON 

/SAVl/ 

BR, BI 




6500 FORMAT! • INTEGRATION DIO NOT CONVERGE - REAL DIFFERENCE 

1 EL3.5,/* IMAGINARY DIFF^^RENCE =* t E13.5 ,5X* THETA =*»E13-5/ 

2 /' Z =*,E13.5,5X'BTR =• , E13. 5 , 5X* BTI =*,E13.5I 

155 FORHATIIH , 5X* 13 , 5X, FIO. A, 2X , • • ,2 X rFLO, 4 ,2 X , • I * ,8 X, • Z = *,F10.4) 

calculation of the 61 MATRIX ELEMENTS 

FMW = FLOAT(MW) 

DO 300 JC=l,N 
IF(JC-I) 310,310,311 

310 ZZ=-C(JC)^B2 
GO TO 312 

311 JF=JC-l 

ZZ=ZZ»(C(JC)+C(JFn*BZ 

312 ICN=ICN+1 
ZZSQ = ZZ * Z7 
HRITE{6,160I tCNfZZ 

160 FORMAT! 1H1,5X, 'BLOCK ND»,I3,5X,*Z DISTANCE = *,F10.4,//» 

IDN=0 

DO 330 K = l,M 
B2D = B2 * OIK) 

IFIX- 1)321, 321,322 

321 R= B2D -t- AINK 

GO TO 323 

322 KL=K-l 

R=^-R + ! D(KL )+OlK ))^B2 

323 IDN=inN+l 

BR! ICN, IDN 1=0-0 
BI ( ICN, IDN )=0.0 
DO 370 I = l,L 




IF 

1 1 - 

2) 

334 

» 

335, 

334 

DT 


DTI 






TT 


0.0 










GO 

TO 

336 

335 

DT 


DT2 






TT 

= 

! DTI 

DT) * .5 









GO 

TO 

336 

337 

DT 

= 

DT2 






TT 

= 

TT 

DT 





336 BDTD4 = B ★ DT / 4,0 
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DTD2 = OT / 2.0 
HOrOAD = BDTDA ♦ DIKI 
DO 350 J=1,LP 
ILK=J 

ILL=NN + 2*( J-1) 

BTR( J)=0-0 
RT I( J )=0.Q 

CALL GAUSSi X,WT, ILL) 

SANR=0.0 

SANI=0.0 

nn 355 I 1=1, ILL 

RI =R- X( I n *B 20 

RISQ = RI ♦ Ri 

DO 355 J J=l, ILL 

TI=TT- X{JJ)*DT02 

Rn = {(B-RI<' C0S(Tn)^*2 + RISQ ♦ ( Sn^(Tin*+2 + ZZSQ 

APT=RnX'HDS 

CART = COSIART) 

SART = SIN(ART) 

SINT=R I*( B-RI* CnSITD) * COS( FMW*TI J /(Ra^RO) 

DUMY = SINT ’J' WTUI) ♦ WT(JJ) 

SANR=DUMY^(CATT/Rn+ SART+WDS) +SANR 

5ANI=0UMY* { SART/RO- CART*WOS) tSANI 

355 CONTINUE 

BTRl J )=SANR*ROTDAD 
BT I( J )=SAN I + RDTDAD 

IF IJ ,LT. 2) GO TO 350 

JL=J“1 

ABT = ABSI ArtSIBTRIJ}) - ABSIRTRUU)) 

IF (ART .GT. CHK) GO TO 357 

ARC = ABS( ABS(BTl(J)l - ABSIRTKJLm 

IF (ABC .LE. CHK) GO TO 360 

357 IF(J-LP)350, 353, 360 

353 WRITE 16,6500) ART, ABC, TT , R ,BTR(J), QTI(J) 

350 CONTINUE 
360 KIL=ILK“1 

BR( ICN, ION )=BR ( I CN , I ON ) +BTR ( ILK) 

BH ICN, ION )=BI( ICNrIDNM-BTIIILK) 

370 CONTINUE 

IF(MnO( IDN,45) .EQ.n) WR I TF ( 6 , 6 843 ) 

6R43 FORMATt IHU 

WRITE! 6, 16 5) I ON , 5R ( ICN, ION) ,R I ( ICN, IDN) , R 
165 FORHATdH , 5X, 13 , 5X, F 10. 4 , 2X , • , 2 X,E1 0, 4 , 2 X, * I • ,8 X , » R = *,F10.4) 

330 CONTINUE 

CALCULATION OF THE Cll MATRIX ELEMENTS 

00 430 K = 1,N 
R2C = R2 ♦ C(<) 
rF(K-l)42l,42l,422 
421 Z=-B2C 
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GO TO A3 3 
A22 KL=K-l 

Z = Z-(C(KLM-C(K ) )*B2 
A3 3 IDN=IDN+l 

ZB = ABS (Z) 

BRI ICN, IDNI=0.0 
BM ICN, IDN)=0.0 
DO 470 1=1, L 

ADDITION TO TEST CLOSED FORM SOLN. 

IF II.GT.l) GO TO 4446 

ILK = 2 

BT U ILK ) = O.DO 
BTR( UK) = O.DO 

IF {K .NE. vtC) GO TO 460 

BTRl ILK) = DTI 

IF (MW ,NE. 0) BTRdLK) = Z-0 ♦ SIN(FMW*DT1 /2.D) / FHH 

GO TO 460 


4446 CONTINUE 

C — - TD HERE - - 

IF ( I - 2) 435 , 435, 437 

435 OT = DT2 

TT = { Dll + DT ) * .5 

GO TO 436 

437 DT = DT2 

TT = TT + DT 

436 8DTD4 = B ♦ DT / 4.0 
DT02 = DT / 2.0 
BDTD4C = BDTD4 ♦ CIK) 

DO 450 J=1,LP 
ILL=NN+2*( J-1) 

ILK=J 

BTRI J )=0.0 
BTI(J)=0.0 

CALL GAUSS! X»WT, ILL) 

SANR=0 .0 

SANI=0.0 

DO 455 I 1=1, ILL 

ZI=Z- X( II ) ♦B2C 

ZZZISQ = 

DABSZI = ABS (ZI) 

DO 455 JJ = 1, ILL 
TI=TT- X(JJ)*0TD2 

Rl = (l cost TI) )^*2 (B* SIN(TI))^^2 ZZZISQ )*^0.5 


ART=RI*W0S 
CART = COSIART) 

SART = SIN(ART) 

SINT=B^(1.0- CDS(TI)) ♦ COS(FHW^TI) / IRI^RI) 

OUMY = R * SINT 4 WTIII) * WTUJ) 

SANR = DUMY*(CARTyRH- SART*MDS) ^-SANR 

S4NI=DUMYt( SART/RI" GART^MOS) 4-SANl 

455 CONTINUE 
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BTR( J »=SANR=*'BDTD4C 
8T M J >=SAN I*BDT04C 

IF (J .LT. 2t GD TD 450 

JL=J-1 

ABT = ABS( ABS(BTR.(jn- AB5[8TRCJL>)> 

IF (ABT -GT. CHK) GO TO 457 

ABC = ABSt ABSIBTKJ)) - ADSIBTKJLMI 

IF (ABC -LE. GHK I GO FQ 460 

^57 IFlJ-LP)450f453,460 

453 rtRlTF (6r6500) ABT, ABC, TT, Z ,BTR(J), BTI(J) 

450 CONTINUE 
460 KIL=ILK-1 

BR{ ICN,IDN)=0R( ICN, lDN»i-BTR(ILK> 

BIC ICN, IDN ) = BI( ICN, lONI^-BTIt ILK) 

470 CONTINUE 

iF(MOO{ inN,45) .EU.O) WR I TF ( 6 , 6843 ) 

WRITFI 6, 166) inN,BR{ ICN, IDN) ,81 (ICNflDN) , Z 
166 FORHATflM , 5X, 13 ,5 X , FIO . 4, 2X , ' <■ * , 2 X ,F 10. 4 , 2X, * I ' , 8 X , * Z = • ,F10.4) 
430 CONTINUE 

CALCULATION OF THE CI2 MATRIX ELEMENTS 

IF( INK-l ) 300,600,600 
600 DO 630 K=l,N2 

B2C = R2 *C2(K ) 

IFlK-l)63l,631,632 

631 7= B2C 

GO TO 633 

632 KL=K-1 

Z=Z+(C2(KLI *C2(K) ) ♦ B2 

633 IDN=IDN+l 

ZB = ABS (Z) 

BR( ICN, ION 1=0.0 
BI( ICN,IDN)=0.0 
DO 670 i=l,L 




IF 

{ t - 

2) 

634 

t 1 

635 , 

634 

DT 

= 

DTI 






TT 

= 

0,0 










GO 

TD 

636 

635 

DT 


0T2 






TT 

= 

( DTL 

DTi ♦ .5 









GO 

TO 

636 

637 

DT 

- 

0T2 






TT 


TT 

DT 




636 

BDTDA 

= R 

♦ DT / 4.0 





DTD2 

= DT 

/ 2.0 





BDTD4C = 8DTD4 <^C2{K) 

DO 650 J=I,LP 

ILK=J 

ILL=NNf2*l J-1) 

8TR( J )=0.0 
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DT I( J )=0.0 

CALL GAUSS( XrWT, ILL) 

SANR=0,0 
SAM 1=0.0 
DO 655 I 1=1, ILL 

zr=z- x{in+B2c 

ZZZISQ = 

DABSZI = ABS (Zn 
DO 655 JJ=1,ILL 
TI=TT- XIJJJ^OTOZ 
CIL = CDSdU 
StL = SINITI) 

RI = [( B-A*CIL ZZZISQ )**0,5 

AR T=Rl^W0S 
CART = cost ART) 

SART = SINIART) 

SINT=l B-A*CIL ) <^A^ COSIFMW^TI) / IRI*RI) 

DUMY = SIMT ♦ WTIII) ♦ WT(JJ) 

SANR=DUMY»ICART/RI+ SART*WOS» 

SANI=0L‘MY*( SART/RI- CART^ WDS) 

655 CONTINUE 

BTR( J)=SANR^BDTD4C 
BT IIJ ) = SAN I+BDTD4C 

IF {J .LT. 2) GO TO 650 


/ 


+/SANR 
+ SANI 


JL = J-1 

ABT = ABS; ABSIBTRU) ) 
IF (ABT ,GT. CHK) 

ABC = ABSI ARSIBTIIJ)! 
IF (ABC .LE, CHK) 

657 IFIJ-LP) 650,653, 660 
653 WRITE 16,6500) ABT, AB 
650 CONTINUE 


ABSIBTR ( JL) ) ) 



GO 

TO 

657 

ABSIBTK JL) ) ) 



GO 

TO 

660 


, TT, Z ,BTRiJ), BTUJ) 


660 KIL=ILK-l 

BR( ICN, IDN ) = BR( ICN , I ON ) +BTR ( I LK ) 

Bl( ICN, lDN) = BinCN, IDNJ + BTHILK) 

670 CONTINUE 

IFIMODl IDN, 45) .FQ.O) WRI TE (6, 6 843) 

WRITE! 6, 15 5) ION,BR(TCN, ION) ,BI(ICN,IDN) ,Z 
630 CONTINUE 
300 CONTINUE 
RETURN 
END 
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SUBROUTINE CAL 700 


C ZDM915 64 - 61 

DIMENSION C1401tDI40) , BR ( 60, 60 1 , BTR ( 60» , B K 60 1 60 ) ♦ BT I ( 60 ) 


1 

fX( 16), 

WT(16) ,C2(40» 



COMMON 

/GAUS/ 




1 

BSQ 

, WDS, ZL, RU, 

ASQ ,DT1, 

0T2 

2 


, ICN 



COMMON 

/MATGUS/ 

B2, 

AINK, PI 2, 

C, D, 

1 » 

INK, N, M, 

MW, LP, Lf NN 

,N2 


COMMON 

/SAVl/ 

BR,BI 




6500 format (* INTEGRATION DID NOT CONVERGE - REAL DIFFERENCE =» r 

1 E13,5, 5X, * IMAGINARY DIFFERENCE =• ,E13.5/15X, *THETA =*, 

2 E13.5, 5X, 'Z , E13.5 , 5X,*BTR =»,FIT.5,5X, • BT I =*,E13.5'J 

155 FDRMATIIH , 5X, 13 , 5X , FIO. 4, ZX , * <- • , 2 X,F 10. 4 ,2 X , ■ I * , 8 X , » Z = *,F10.4) 
160 FDRMAT(1H0,5X, 'BLOCK ND',I3,5X,»Z DISTANCE = *,F13.4,//) 

165 FOkTATIIH ,5X, I3,5X,F10.4,?X,*4-' ,2X,F10.4,2X,' I* ,8X,*R = »,F10.4> 

166 FORMATIIH , 5X, 1 3 , 5 X, F 10. 4, 2X , • ♦ ' , 2 X,F 10 . 4 , 2X , ' I • , 8 X , ' Z = ',F10.4} 
C + 

FMW = FLOATCMW) 

DO 700 JC=1,N2 
IF{JC-11710,71Q, 711 

710 77=C2(JCI * B? 

GO TO 712 

711 JE=JC-l 

ZZ=ZZ4-{C2IJC) + CZtJFl) * B2 

712 [CN=ICN+1 
ZZSQ = 17 * ll 
WRITFI6,160) ICNfZZ 
IDN-0 


CALCULATION OF THE B2 MATRIX ELEMENTS 


DO 730 K=1,M 
B2D = B2 * D(< ) 
IFIK-D721, 721,722 

721 R = B20 + AINK 

GO TO 723 

722 KL=K-1 

R=RttO(KI + n(KL n*B2 

723 IDN=ION-H 

BR( ICN, IDN)=0 .'> 

BH ICN, IDN »=n. V’ 

DO 770 I = l,L 



IF 

(1-2) 

734, 735, 

734 

DT = 

DTI 



TT = 

0.0 

GO TO 736 

735 

DT = 

DT2 



TT = 

(DTI + DT) ♦ ,5 

GO TO 736 

737 

DT = 

DT2 
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TT = TT + DT 

736 BDTD4 = B ♦ DT / 4.0 
DTD? = DT / 2. 0 
RDTD4D = BDTD4 ♦ DtK) 
on 750 J=1»LP 
ILK=J 

ILL=NN<“2*( J-ll 
BTR< J >=0.0 
BTK J > = 0.0 

CALL GAUSSrXtWT, ILL) 
SANR=0.0 


SANI=0.n 
DO 755 I 1=1, ILL 
RI=R- X{ II )*B2D 

RISQ = RI * RI 
DO 755 JJ=1,ILL 
TI=TT- XIJJ)+DTD2 

RQ = ((A-RI* CQS(TI>)*=«'2 * RISQ *( SIN(TI)>**2 + ZZSQ >♦♦0.5 


ART=RQ*WDS 
CART = COSIART) 

SART = SIN(ART) 

SINT=RI*(A-RI* COS(TI)) * COS(FMW*Tl) / (RC*RO> 
DUMY = SIMT ♦ WT(in ♦ WTIJJ> 

SANR=0UMy<'{CART/R04- SART*WOS> +SANR 

5ANI = DUHY»I SART/RO’ CART+HDS) ♦•SANI 

755 CONTINUE 

BTR( J l=SANR*0DTD4n 
BT H J >=SANl*BDT0 4D 

IF (J .LT, 2) GO TO 750 

JL=J-1 


ABT = ABSt ABSIBTRUn - ABS(BTRIJl))) 

IF (ABT .GT. CHK> GG TO 757 

ABC = ABSI ABStSTKJl) - ABSIBTKJLH) 

IF (ABC .LE. CHK> GO TO 760 

757 IF(J-LP)750,753, 760 

753 WRITF (6,6500) ABT, ABC, TT, R ,BTR{J>, BTIIJ) 
750 CONTINUE 


760 KIL = Il K-1 

BR( ICN, IDN )=BR( ICN,IDN>+BTR{ ILK) 

BIl ICN, IDN )=SI ( ICN,I0N> + BTI I ILK) 

770 CONTINUE 

WRITE( 6, 16 5) ION ,BR( ICN, IDN) ,BI ( ICN, IDN) ,R 
730 CONTINUE 


CALCULATION OF THE C2l MATRIX ELEMENTS 

DO 830 K=1,N 
B2C = B2 ♦ C(K> 

IF(K-1)821, 821,822 
821 Z=-82C 

GO TO 833 
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822 KL=K“1 

Z;=Z-tC(Kl+C(KL ) )*B2 
833 I!1N=IDN + 1 

78 = ABS (Z) 

ZZSQ = 11 * ll 
8R( ICM, IDN )=0.0 
QIl ICN, IDN )=0.0 
DO 870 1=1, L 




IF 

I I - 

2) 


834, 835, 

834 

DT 

= 

OTl 





TT 

— 

0.0 



GO TO 836 

835 

DT 


DT2 





TT 

- 

1 DTI 

+ DT » * .5 


GO TO 836 

837 

DT 

= 

DT2 





TT 

■= 

TT + 

DT 




836 BDTD4 = R ♦ OT / 4.0 
DTD2 = DT / 2.0 
BDTD4C = B0T04 * C(K) 

DO 850 J=1,LP 
ILL=NN+2*l J-1) 

ILK = J 

BTR( J }=n,o 
RTIt J »=0.0 

CALL GAUSS{ X,WT, ILL) 

SANR=0,0 
SANI=O,0 
DO 855 1 1=1, ILL 
ZI=Z- X( II )*82C 

77ZISQ = (ZZ-Zn**2 
DABSZI = ABS (ZI) 

DO 855 JJ=l,ILL 
TI=TT- XtJJ)*0TD2 

RI=(IA-B* CClSm))**2 + BSQ ♦ ( SIN(TIU**2 *■ ZZZISO 1**0.5 

ART=RI^W DS 

CART = CO St ART) 

SART = SIN(ART) 

SINT = R*tA-B* COSITin * COSIFMH^TI) / (RI*RI) 

DUMY = SINT ♦ HT(Il) ♦ WT(JJ) 

SANR = DUMY<‘(CART/RH- SART+WDS) +SANR 

SANI = 01JMY*^I SART/Rl- CART»WDS1 +SANI 

855 CONTINUE 

BTRt J) =SANR*BDTD4C 
13TK J )=SANI»BDTD4C 

IF tJ .LT. 2) GO TO 850 

JL=J-1 

ABT = ARSt ABStBTR(JJ) - ARSCBTRULin 

IF (ABT .GT. CHK) GO TO 857 

ABC = ABSl ARSlBTHJn - ABSIRTIULD) 

IF (ABC .LE. CHK I GO TO 860 
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857 

853 

850 

860 


870 


830 


BTRIJJ t BTKJ) 


C 

c 

c 


921 

922 
933 




IF( J-LPI850,853, 060 
WRITH (6,6500) AST, ABC, TT, Z 
CONTINUE 
KIL=ILK-1 

BRt ICN, IDN )=BR{ ICN,IDN)4-BTR( ILK) 

BII ICN, lONJ^BIt ICNflONH-BTHILK) 

CONTINUE 

WRITE! 6, 166) IDN,8R! ICN, ION) ,8 I (ICN, ION) ,Z 
CONTINUE 

CALCULATION OF THE C22 MATRIX ELEMENTS 

DO 930 K=l,N2 
B2C = B2 *C2(K) 

IFIK-U921, 921,922 
Z = B2C 
GO TO 933 
KL=K-1 

Z=Z + (C2(K) *■ C2(KD) ♦ B2 

IUN=IDN+l 

ZB = ABS (Z) 

BR( ICN, I0N) = 0.0 
8U ICN,IDN )=0.0 
00,970 1=1, L 

ADDITION TO TEST CLOSED FORM SOLN. 
IF ( I.GT.l) 

ILK = 2 

8THILK) = 0,00 
BTR(ILK) = O-DO 
IF (K .NE. JC) 

BTR(ILK) =-DTl 


GO TO 4446 


GO TO 960 



IF 

(MW . 

.NF. 0) BTRIILK) 

=-2.0 * SIN(FMW*DTl/2.0» / FHW 





GO TO 960 

4446 

CONTINUE 




— 

- 

TO HERE - 

----- - - - 


IF 

( I - 

2) 

935, 935, 937 

935 

DT = 

DT2 




TT = 

(DTI 

3- DT) * .5 






GO TO 936 

937 

DT = 

0T2 




TT = 

TT 

DT 


936 

BDTD4 = B 

♦ DT / 4,0 



DT 02 

= DT 

/ 2.0 



8DTD4C = BDT04 *C2(K) 

DO 950 J=l,LP 

ILK=J 

ILt=NN«-2*( J-1) 

BTR( JI-O.O 
BT I( J)=0.0 

CALL GAUSS! X,WT, ILL! 
SANR=0*0 


S 111 wip» 
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i 


SANI=0,0 


DO 955 I 1=1, ILL 
ZI=7- X( in* 82 C 

ZZZISQ = [ZZ-ll\**Z 
DABSZI = ABS (Zn 
DO 955 J J=1 , ILL 
TI=TT- XIJJJ*DTD2 

RI=UA-A* CQS(Tin**2 + (A ♦ SIN(TI))**2 + ZZZISQ 

SINT=A*(1.0- COSITin ♦ COS(FHW<'TI) / (RI*RI) 

ART = Rl ♦ WDS 
CART = CnSIARTI 
SART = SIN(ART) 

DUMY = A * SHT + WT(IT) * WTtJJ) 

SANR=OUMY* ICART/RH- SART^WDSI +SANR 

SANI=OUMY* (5ART/RI- CART*WDS) +SANT 

955 CONTINUE 

BTRtJ J=SANR*BDTDAC 
BTIU )=SAN !*R0TD4C 

IF (J .LT. 2) GO TO 950 


Jl =J-1 


AHT = ABSI ABS(t3TR(JM - ABSCBTRULlU 
IF lABT .GT, CHK) 

ABC = ABSI ABSIBTUjn - ABSlSTItJUM 
IF (ABC •LE, CHKJ 
957 IF( J-LP)950,953,960 

953 WRITE (6,6500i ABT, ABC, TT, Z ,8TR{J), 
950 CONTINUF 
960 K1L=ILK-1 

BR ( !CN, IDN )=RRI ICN , I DN HRTRI ILK I 
BI ( ICN,10N) = BI{ ICN,IDN»*-BTI( ILK) 


GO TO 957 
GO TO 960 
BTK J) 


970 continue 

WRITE! 6, 15 5) IDN,RR( ICN, I DN» f B I IIC N, IDNI ,Z 
930 CONTINUE 
700 CONTINUF 
RETURN 
END 
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) **0 . 5 




APPENDIX 4 


GAUSS SUBROUTINE LISTING 






ZDM915 11 - BT 


SUBROUTINE GAUSS ( X,WT, M) 

C STORE DP TEST IN 81 

CD IMPLICIT REALMS CA“H,D~Z) 

DIMENSION X( 12),WT(121 
IFIM-4H0, 11,12 

10 X( U=-0. 677350269189626 
X( 2)= -XC 1) 

WTU) = 1.0 
WT(2)=l,0 
GO TO 20 

11 X( lJ=-0, 861136311594053 

X{ 2I=“0,339981 043584856 
Xt3)=“ X12) 

x(4j=- xm 

WT( 1 1=0.347854845137454 
WT(2)=0.65214S154862546 
WT(3)=WT(2* 

WTt4)=HT{ n 
GO TO 20 

12 IFIM-8U3, 14,15 

13 XI lI=-0. 932469514203152 

X( 2)=-0.661209386466265 
X( 3l=-0. 238619186083197 
X(4I=- X(3) 

X(5)=- XI2) 

x{6»=- xm 

vJTIl ) = 0. 1713244923 79170 
WT(2»=0. 360761573048139 
WTI 31=0.467913934572691 
WT(4t=WT(3) 

WT(5I=WT( 21 
WTl6)=WTm 
GO TO 20 

14 XI 1)=“0. 960289856497536 

XI 2)=-0. 796666477413627 
XI 3I=-0, 52 5532409916 329 
XI 4)=-0, 183434642495650 
XI5»=- XE4) 

XI6J=- X(3) 

X(7I=- XI2> 

X(8)=- XIU 

WTI l) = 0, 101228536290376 

WTI 21=0.222381034453374 

'*<T( 31 = 0,313706645877887 

WTI 41 =0.362683 783378362 

WTI 5>=WT14J 

WTI61=WTI 3) 

wn7) = WTI21 

WTI8)=WT( 1) 

GO TO 20 

15 IFIM-121 16, 17, 18 

16 XI l)=-0. 97390652851717? 

XI 21=-0. 865063366688985 
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X( 3}=-0. 67940956 8Z99024 
X( 4) =-0.433395394129247 
X{ 5)=-0. 14B874338981631 
X(6»=- X(5) 

X{7>=- X(4) 

X18>=- X13I 

X(91=- X(2t 

X(10)=- X(l) 

WTU» = 0.C6667l 344308668 
WT(2»=0. 149451349150581 
MT{ 3 1=0,219086362515982 
WT {41 = 0.2692 66 71 93 09996 
vITt 5) =0.295524224714 753 
WTt6)=WT(51 
WT17>=WT(4» 

WT( B)=WT( 3) 

WTI9)=WT(21 
WT(10»=WT( 1) 
on TO 20 

17 X{ 1)=-0.981560634246719 

X( 2>=-0. 904117256370475 
X( 3)=-0, 769902674194305 
X( 4)=-0. 587317954286617 
XI 5)=-0. 367831498918180 
X( 6»=-0. 125333408511469 
X171=- X(6J 

Xt81=- X(5) 

X(9)=- X{41 

X(10>=- X(3J 

X{ IH=- X(2) 

xti2)=- xm 

WTt l)=0. 047175336336512 
WT( 2) =0.106939325995318 
WT(3»=0. 160078328543346 
WT( 4) =0.203167426723066 
WT(5»=0. 233492536538355 
4T(6)=0, 249147045813403 
WTI7)=WT16) 

WT( 0J=WT{5) 

WT(9)=WT(4) 

WTl 10)=WT( 31 
WTU1I=WT( 2) 

WTt 12)=WT( 1) 

GO TO 20 

18 IF lM-16) 19t 2lf 20 

19 Xt ll=-0. 986283808696812 
XI 21=-0, 928434883663574 
X( 31=-0. 827201315069 765 
XI 4l=-0. 687292 904811685 
X[ 5)=-0. 515248636358154 
XI 6)=-0. 319112368927890 
X( 7)=-0. 108054948707344 
X(8)=- X(7t 
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X( 9|=- 


XI 6) 

X( 101 = 


-XI 5) 

X( 11) = 

- 

XI4) 

X( 121 = 

- 

XI3) 

X( 131 = 

- 

XI2 ) 

XI 141 = 

- 

XII ) 

WT( 11 = 

0 

.035119460331752 

WT121= 

0 

.080158087159760 

WT{3)= 

0 

. 121518570687903 

HT(4)= 

0 

.157203167158194 

WT15)= 

0 

. 185538397477938 

WT(61 = 

0 

.205198463721296 

WT(7)= 

0 

.215263853463158 

WTl 81 = 

WTl 7) 

WTI91 = 

WTl 6) 

WTl 10) 

- 

WTl 5) 

Wn 11) 

- 

WTl 4) 

WT(12I 

= 

WTI31 

WTl 13) 

= 

WTl 2) 

WTl 14) 


WTl 1) 


21 XI 1) =-0.9894009^4991650 
X( 2)=-n. 944575 02 30 73233 
X{ 3 J =-0.865631 202387832 
X! 4>=-0. 755404408355003 
X( 5»=-0, 617876244402644 
XI 6)=-0. 458016777657227 
XI 7)=-0. 281603550779259 
X{ 8»=-0. 09 5012509837637 
X(9»=- X(8) 

X110J= -X(7) 

XtllJ=- X{6) 

X{ 121 = - X(5) 

XC13)=- X(4I 

X(14)=- X(3) 

X(15I=- X(2) 

XM61=- XdJ 

WT{ 1) = 0. 027152459411754 
WT( 2 1=0.062253 523938648 
HT( 31=0.095158511682493 
WT( 4) =0.124628 97 1255534 
WT(5)=0. 149595988816577 
WT(6I=0. 169156519395003 
WTl 7 1=0. 182603415044924 
WT(8)=0, 189450610455069 
WTC91=WTI 81 
WT(ini=WT( 7) 
wTmi=wn 6) 

WTt 12)=WTI 51 
WT(i31=WT(41 
HT(14»=WT( 31 
WT( 151=WT( 2) 

WTl 161=WT( 1) 

20 RETURN 
END 


GO TP 20 
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MATX SUBROUTINE LISTING 
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SUBROUTINE MATX CDRECTP, PUNCH, SAVTP, TITLE, DUMEE) 

C MATRIX CALC FROM 

C BOX METHOD INTFGRATIQN FOR SOURCE DISTRIBUTION OVER DUCT 

C DOUBLE PERC. TEST ZDM915 66 - 66 

implicit RFAL*8 (A-H,0-7» 

LOGICAL ♦ 1 PUNCH 
INTEGER ♦ 2 DATE! 3J 

INTEGER *4 DRECTP , SAVTP ,T!TLE(18J , TPS/07/ 

1 ,MXM/30/, HXN/40/, MXN2/30/, MXNN2/45/ ,MXNUM/60/ 

2 ,RAD{3l/'4*» RADIAL «/, I DUC T ( 31 /• I NNER DUCT •/ 

3 fODUCTt 3)/*OUTFR DUCT »/ 

REAL *4 AS, BS, WS, CHKS, DSIU, CS(U, ESIl), PHIISIU 

1 tPHizsm, vRsm ,ss ,c 2 sm 

COMPLEX *16 VI40J 

1 ,PHI l,PHI2*AA,BB,CC,DD,DCMPLX ,DZERO/I 0. DO, 0. 00 ) / ,DCONJG 
DIMENSION CI40),DU01 ,C2I401 

2 ,PHI1(50) ,PH12(40» ,VSI1) 

3 ,AA{40,45),BB(46,40) ,CC (60,60) ,DD(60,60I 
COMMON /MATRIX/ 

1 V, PHIl, PHI2 ,W, S, E(40) 

2 , IV ,KDF 

COMMON /MATGUS/ B2, ASNK, PI2 ,C ,D ,A ,B ,CHK ,C2 

I ,INK, N, M, MW, LP, L, NN ,N2 

COMMON /SAVE2/ BB 

COMMON /SAVF3/ AA 

C - - CC AND OD MUST BE IN THIS ORDER TO WORK WITH INVT ♦♦***♦ 

COMMON /SAVl/ CC,DD 
EQUIVALENCE I V 1 1 ) , VS 1 1 ) , VRSI 11 ) 

1 ,(A,AS) ,( B,BS1, (W,WS) ,(S,SS), <CHK,CHKS) ,(Om,DS(L)) 

2 ,( C( 1) ,CS( 1) ) ,{Efl),ESm), IPHIini, PHIISIU 1 

3 ,(PHI2t llfPHIPSi U3 ,(C2(11 ,C2Sm» 

100 F0RMATI3612) 

102 F0RMAT(7F10-3) 

103 FORMAT! 18A41 

104 FORHATITIlOl 

105 FURMAT(6E12,5) 

106 FORMAT 1*1', lOX, 18A4//1 

107 FORMAT UU, 20X , 'CALCULATE 

I 7X, A2, 2!*/*,A2J, 5X, 

108 FORMAT t *0 C BUFFER (-Z1*/ 
no FORMAT I'O C2 BUFFER ( + Z1V 
loq FORMAT (*0 D BUFFER !Rl '/ 

200 FORMAT !6E 12-51 

9007 FORMAT I1P6D12.5) 

9008 FORMAT !A2, 2!*/*,A2), 18A41 
C* - 


PHU'S FROM BOX INTEGRATION MATRIX*/ 

18A4///1 

( 8X,1P4D15.4) 1 

iaX,lP4D15.4l I 

(8X,lP4D15-4) 1 




- * - 


CALL DATEIS (DATE! 

WRITE !6,107l DATE, TITLE 
REWIND SAVTP 

RFAD tSAVTPI INK, N, M, MW, LP, L, NN, KDF ,N2 
1 ,B2, AINK, PI2, A, B, CHK, U, S 
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2 iCCU)fi = lTNJ , (OdUI^lfH) tCHnWI^ltU , (C2 (I I , I- UNZ I 

WRITFt6,5189J INK, N, M, MW, LP, L, NN, KDF, N2 , B2, A^NK, 

1 PI2, A, B, CHK, U, S 

5189 FORMAT!* INK=*,I3,« N=*,I3,» M=',I3,* HW=*,I3,* LP=*,I3, 

1 * L=»,I3,* NN=*,I3,» KDF-»,I3,» N2=*,I3/» B?=*,F6.3, 

2 * AINK=*, F6.3, * PI2=*,F9.5,* A=*,F6.3,» B=*,F6.3,' CHK=* , 

3 F8.6/* W = ', FIO.4, » S=* , F10.4I 

WRITE I6,108» IC(U,I = 1,NI 

IF I INK .GT. 0) WRITE (6,110» IC2 tl I , I -1 ,N2I 

WRITE I6,109» !DI I), 1=1, Ml 
IF !.NDT. PUNCH! 


IF I INK .GT, 01 

1003 NO = N 

NN2 = NO + N2 
NUM = M + N + N2 
IF {M.GT.MXM .OR. 

1 NN2.GT.MXNN2 

N02 = 2 * NO 
N04 = 4 ♦ NO 
N22 = 2 * N2 
N24 = 4 * N2 
HR2 = 2 ♦ M 
MR4 = 4 ♦ M 
LR2 = 2 L 
IF! INK-1 )87, 77,77 
77 N = N + N2 
87 CONTINUE 

READ !SAVTP) IIDOIl.Jl, I=l,NUMI, J=1 ,NUM! 

DO 10 1=1, M 
DO 10 J=1,N 
K=J+H 

AA! I,J I = 00! I,K> 
in CONTINUE 
DO 12 I=l,N 
DO 12 J=1,M 
K=I+M 

BBI I,J ) = DDIK.JI 
12 CONTINUE 
DO 15 1=1, N 
DO 15 J=1,N 
K= I+M 
L=J+M 

CCIItJ I = DO!K,L I * P12 
15 CONTINUE 
DO 50 1=1, N 
DO 50 J=1,N 
DD! I, J I = DZERO 
DO 55 K=l,M 


WRITE ITPS,9008) 
WRITE (TPS, 9007) 
WRITE (TPS, 90071 
WRITE (TPS, 90071 


UAlb, UILt 

(cni ,I=L,NI 
(C2!I) ,l=l,N2) 
(0(11, 1=1, HI 


N.GT.HXN .OR. N2.GT.MXN2 .OR. 

.□R. NUH.GT.MXNUM) GO TO 9080 






-v.‘. v; 
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55 DOC I,J > = BB( IfK »*AA(Kt J)*ODU ,J1 
50 CONTINUE 

CN = n.n 
DO 56 I = 1,N 
DO 56 J=l,N 

DOt I,J J = CC( I,J i-DDII ,JI 

CN = CN + DDU,JJ * DCONJGCDDCIfJ) I 

56 CONTINUE 

C — — OD IS IN COMMON AND IS DESTROYED BY INVT *♦*♦♦**♦ 

CALL INVTtDD.DCfNtKKKi 

IF (KKK .E3. U GO TO 9000 

CM = 0.0 
DO 60 I = ltN 
00 60 J=1,M 

CM = CM + CCtl.J) ♦ OCONJGCCCC I# JU 
DD(I,J) = DZERO 
DO 66 K=1,N 

ODdfJ ) = CC( I,K )*BQCK,JKnD(I, J) 

66 CONTINUE 
60 CONTINUE 

COND = DSQRT { CN * CM) / DFLOAT(N) 

WRITE t 6, 7001) COND 

7001 FORMAT (=0»,20X, ’CONDITION NO. =• , 1PE17.8) 

C - - - OUTPUT ALL INFO TO THE DIRECTIVITY FACTOR AND 

C - - - IMPEDANCE PROGRAMS IN SINGLE PERCISION - 

IF (PUNCH) WRITE (7,100) K, M, N, NZ , INK, LP, MW, NN, L 

IF (PUNCH) WRITE (7,200) A, B, W, S, CHK 

IF (PUNCH) WRITE (7,200) (DSC I ) , I =l , HR2 ,2 ) 

IF (PUNCH) WRITE (7,200) ( C SII ) , I = 1 , N02 ,2 ) 

IF (PUNCH .AND. INK.GT.O) WRITE (7,200) (C2SI I ) , 1=1 ,N22 ,2) 

IF (PUNCH) WRITE (7,200) ( E S( I ) , I =1 , N02 , U 

READ IN VELOCITIES 
READ(5,2000,EN0=9000) IV 
2000 FORMAT(2I2) 

C - - IV .GT. 0 OLD CARD INPUT 

C - - IV -E3. n CALL VELINT 

C - - IV .LT. n READ CHAMBER’S OUTPUT 

IVF = IV 
NAP == 1 
NAM = I 

IF (IV) 1005, 1008, 1010 

1005 READ (5, 103) TITLE 
WRITE (6,106) TITLE 
READ (5,104) NAP, NAM, NRM, NSR 
WRITE(6, 104) NAP, NAM, NRM, NSR 
IV = NRM 

GO TO lom 

1008 CALL VELINT (IV, M, D, INK, A, B, AINK, MW) 

1010 M2 = M * 2 

DD 2001 INAP =1,NAP 
DU 2001 INAM =1,NAH 



y 

no 



8 ? 


IF{ 1-1)82, 82,83 
7=- C( I ) * B2 
GO TO 181 
83 KI=I-1 

Z =7 - 1C (KI) * C tl ) ) *82 
181 CONTIMUE 
KOK=I+M 

IF(«QO(KOK,5l) .EQ.O) WRITE (6, 1061 
81 WRITE(6,76) PHlim,Z 

IF (PUNCH) WRITE (7,200) { PHI 1S( I ) ,1=1 ,Na4,2l 

IF( INK-l )97,98,98 

98 00 99 I=l,N2 
IF( 1-1)72, 72,73 

72 Z = C2( I ) * 82 
GO TO 74 

73 Kl=I-l 

Z = Z ♦ (C2IKI) *■ C2(IH ♦ 82 

74 J=I+N(1 

99 CONTINUE 

76 FORMATIIH , 2F10.4,* I*,5X,*Z = »,F10,4I 

WRITE(6,76) PHlllJltZ 
N1 = NO +1 
NR 12 = NQ4 ♦ 1 
NR 22 = N04 + N24 

IF (PUNCH) WRITE (7,200) ( PHI 1 S ( I ) , I =NR12 , NR22 ,2 ) 

97 CONTINUE 
2001 CONTINUE 
9000 RETURN 

9090 WRITE (6,9091) M, NSR 

9091 FORMAT ( *0 THE NO, OF RAOII OD NOT MATCH, CHECK INPUT CAROS* 

1 ,* AND DISK, M I3,9X,*NSR =*, 13) 

GO TO 9000 

9080 WRITE (6,9081) MXM, MXN, MXN2, MXNN2 ,MXNUM 

9081 FORMAT (*0 M IS GT*, 13, *, N IS GT*,T3, *, N2 IS GT*,I3 

1 ,* N+N2 IS GT*, 13, *,nR HtN+N2 IS GT*,I3 / 

2 12X, *THE DATA OP DIMENSIONING MUST* 

3 ,* BE CHANGED* , lOX, 30( **-*)) 

GO TO 9000 

END 
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DO 2001 IRA = li IV 



IF { IVF) 

1015, 1018, 1030 

1015 

IF { INK -EQ- Ot 

GO TO 1016 


IF (NSR -NF- Ml 
NB = 1 

GO TO 9090 
GO TO 1017 

1016 

CONTINUE 


- 

FAN DUCT 

— « - _ - 


IF (NSR -GT. Ml 
NB = M - NSR 

GO TO 9090 


IF {NB -EQ. 01 
00 1116 I =1,NB 

GO TO 1117 

1116 

V( II = 07ER0 


1117 

NB = NB 4- 1 


101? 

READ {5,105) (V(n,I = NB,M) 

GO TO 1040 

1018 

CALL VELCAL lIRA, VS) 
DO 1020 1=1, M 
K = M - I + 1 
J = M2 - I ♦ 2 + 1 
VSUI = VSIK) 


1020 

VSIJ+1) = C.DO 

GO TO 1040 

1030 

00 550 1=1, M 


550 

V{ n = O.DO 



READ(5,102 ,EMD=0000) I VSl f) , I =l t M2 »2 I 

1040 DO 65 1=1, N 

PHIK I I = DZERO 
OQ 70 K=1,M 

70 PHIK H = nD( I,KJ*VtKl»{-l.O)+PHIl(n 
65 CONTINUE 
DO 75 I=l,M 
PHI2{ I I = DZERO 
DO an K=i,N 

80PHl2m=- AA{ I,K)*PHI KK.I/ PI2 +PHI2(I» 

PHlZm=V(I»/ PI2 >PHI2(I» 

75 CONTINUE 

WRITE! 6, 95 » IRA 

95 FORMATCIHI , 'SOURCE STRENGTH 0 1 STP I BUTI ON* , Il5 , /) 

DO 90 1=1, M 

IFlI-1191,91,92 

91 R=DU1* 82 «• AINK 

GO TO 90 

92 KI=I-1 

R=RfrI DI I I + OtKin* 82 
90 HPITE{6,96» PHIZdl.R ,VCII 

96 FORMATIIH , 2F10,4,* I»,5X,*R * ',F10.4, 7X,*V = •♦2F8.2I 

IF {PUNCH* WRITE *7,2001 I PHIZ SI I) ,I = 1 ,HR4,2) 

9009 FORMAT I 3A4, 'PHS'^S FOR ANG. PT»,I3, ANG- HD', 13 

1 , •, AND RAD- MD', 13, / UP6E12.5)) 

DO 81 1 = 1, NO 



c 


SUBROUTINE INVTI A, AI NV,N,KKKJ 


Z0M915 65 - 65 


IMPLICIT REAL*8 (A-H»0-Z) 

COMPLEX *16 A, AINV.S,B,Tf8A0 
COO COMPLEX *08 A, AlNVt S T ,BAO 

DIMENSION M60,60» ,AINV{ 60,601 ,B(60 ,1201 ,BDt2» ,S(120) 

COMMON /SAVl/ B 
EQUIVALENCE (BDUItBAO) 

1000 FORMAT <601,5.41 

C ♦** *♦* ** ♦ 

MM=2*N 
KKK - 0 
00 23 1 = 1, N 
DO 23 J=1,N 

23 B( I,J ) = A( 1,J J 
K = N+1 

DO 24 I=1,N 
on 24 J=K,MM 

24 B( I, J I = <0.0, 0-0) 

DO 25 1=1, N 

K = I + N 

25 B( I,K1 = i 1.0, 0.0) 

DO 33 1=1, N 

J = I 

8AD=B< I, J1 
BR = BD< 1) 

BI = B0(2) 

aAR=[BR*BR+BI*BI 1**0. 5 
IF(BAR.GT.2.0D-06) GO TO 26 
L = I + 1 
C DO 12 K=< ,N 

K = L 

BAD=BIL, J) 

8R = BD( 11 
BI = 8012) 

BAR=< BR*BR + BI*HI |**0.5 
IF(BAR.GT.2.00-06) GO TO 27 
C 12 CONTINUE 
KKK = 1 
WRITE <0,10101 

1010 FORMAT (*0 A SINGULAR MATRIX WAS INPUT TO SUBROUTINE INVT-M 
RETURN 

27 DO 13 H=1,HM 

13 S(M1 = B<K,M) 

DO 14 M=1,MM 

14 BtK,Ml = B(I,M1 
00 16 M=1,MM 

16 B( I, Ml = SIMl 

26 T = B( I, J1 
on 17 K=l,MH 

17 B( l,K) = BII,K l/T 
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ft ^ ni J. 




L = I + 1 

IF IL.GT.NJ G3 TQ 33 
DO 10 K=L,N 
T = 8«K,J) 

DO 10 NN=1,MM 

10 BIK,NN) = B<K,NN)-T*Bn,NN» 
33 CONTINUE 

IF IN.EQ*1I G3 TO 50 
00 20 I-2,N 
J = I 
L = I - 1 
T = B(L,J) 

DO 21 K=J,MM 
21 B(L,K) = B(L,K 

IF (L.EQ.l) GO TO 20 
L = L- 1 
GO TO 29 
20 CONTINUE 
50 DO 30 1=1# N 
DO 30 J=l,N 
NN = J + N 

30 AINVl If J) = BUiNNi 
RETURN 
END 


APPENDIX? 


VELINT SUBROUTINE 


pbiobeing page blank not BnatEO 
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SUBROUTINE VEL INT (NNZ«NARGr RARG, INKt AA, BB, AINKf MW) 

C ZDM915 25 - 25 

- - '^ODAL VELOCITY CALCULATION 

IMPLICIT REAL* BIB I 
DOUBLE PRECISION 
1 ,DABS,DFLOAT 

,BJl,BJ2 fB J,Z,R00TI20) ,C ( 2 01 , 8P,X(At) ) 
,SIG,0N,BY1,BY2,8JS1,BJS2 r8YSl,ByS2fFtFPfZS,XS 
tRARG*R,XKMN.OMN,AA ,AINK 
DIMENSION RARGdl »R(301,BESMm 
1 FnRMATl6E12*6) 

C-* - *-*-*-*-* 

READ(5,3» NNZ 
IFIMM.LT.O) GO TO 11 

3 FORMAT! 3612) 

S IG=AA/BB 
ON=DFLOATI MW) 

C{ U=ON 

IFtMW.EQ.O) Cl 1»=3.D0 

DO 30 J=1,NNZ 

MP=MW+l 

NW=J-1 

N=l 

X( U=C( J.) 

10 CONTINUE 
CZ=UOO 

4 FGRMATdH , 3d OX ,025. 16) 1 
WRITE!6,4) C(JI 
IFIMW.FO.O) GO TO 15 

CALL DBESJYl X( N ) ,MW, B JM, BYM, BJMl ,8 YMl , I FRR ) 

XS=XIN )*SIG 

CALL DBESJY{XS,MW,8JSM,BYSM,B JSM1,8YSMI ,JERR) 
BJPM=-BJH1+DN*BJM/X(N) 

RJPSM=-SIG*BJSHH-DN*BJSM/XIN) 

RYPM=-RYMH-ON*SYM/Xl N) 

8YPSM=-SIG*BYSMU-DN*BYSM/XIN) 

BJPPM=8JM1/X(N 1-lDN/CXlN )*X!N) 1-(DN/X(N) ) **2*1 . DO 1 *B JM 
BYPPM=BYM1/X(N |-{DN/IXIN I* X( N)l - (D N/X ( N) ) **2+ 1 . OO ) *BYM 
BJPPSM=BJSMl*SIG/XlN)~lDN/(XlN)*X!N)l-{ON/XlN) J**2<-SIG**2)*ajSM 
BYPPSM=RYSM1*SIG/X(N)-(DN/IXIN)*XINI )-(DN/X(N) ) **2 +-S I G**2 ) *BYSM 
F= BJPM*BYPSM-BJPSM*BYPM 

FP=BJPPM*0YPSM+BJPH*BYPPSM-8JPPSM*8YPM-BJPSM*BYPPM 
WRITF!6,4» F,FP 
XtN+l)=X{N)-F/FP 
WRITE(6,4) X(N+1I 

CALL DBESJYIXI N*l) ,MW,BJM,BYH,BJH1,BYM1,IERR) 

XS=X(N+1)*SIG 

CALL DBESJY!XS,MW,BJSM,BYSM,BJSM1,BYSMI,IERR) 

BJ=BJPH*BYPSH-BJPSH*BYPM 

IF (N.GT.25) GO TO 20 

IFIOABS! BJ I .LT.l.D-15) GO TO 20 
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GO TO 16 

15 IF(NH.EO,0» GO TO 100 

CALL DRPSJYt XC N J , 1 ,B J 1 ,B Y1 ,B J2 *BY2 ♦ I ERRJ 
WRITE(6,A) BJltBYJ. 

WRITei6,4) BJ2,BY2 
XS^SIG^^XtN ) 

CALL DBFSJYIXS,1,BJS1.BYSI,BJS2,BYS2 tIERR) 

WRITE(6,4) BJSl.BYSl 
WRITE(6,A) BJS2fBYS2 
F^S!G*( BJ1#BYS1-BJS1*BYI J 

FP=SIG*BJ1*(-S IG*RYS2+BYS1 /Xt N»)+SIG*BYSl ♦ I-BJ24-BJI/XC Nll-SIG* 
ABJS1*(-BY2*BY1/X(NI)-SIG*BY1*(-SIG*BJS2<-BJSI/X{N1 I 
X{ N+l)=XtNJ-F/FP 
WRITF(6,4J XlN + n ,F,FP 

CALL DBESJYt X( N + n tl tBJl ,BYl , B J2 ,B Y2 , I ERR J 
XS=SIG*XtN+l) 

CALL OBESJYC XS, 1 ,B JS 1,BYS1 tB JS2 t BY S2 , I ERR I 
BJ = STG*( BJ 1*BYS1-BJS1*BYU 
WRITE{6, 4) BJ 

IF tDABSt BJ » .LT. 1.0-15} GO TO 20 

16 CONTINUE 
N=N+l 

GO TO 10 
100 X( 21 = 0.00 
20 ROOTC J )=X(N+l) 

IFUMW.EQ.O).ANO.(NW.EQ.On GO TO 17 
IF(J .FQ.NNZ } GO TO 30 
IN0N=0 
INDP=n 
1 = 0 

60 1=1+1 

Z = ROOTtJ ) + .2D0*nFL0AT( I ) 

WKITE(6,A) Z 
IFIMW.EQ.O} GO TO 80 

CALL 0BFSJY<Z,MW,BJH,BYM,BJMl,8YM! ,IERR) 

2S=SIG*Z 

CALI OBESJYtZS,MW,BJSM,B YSM,BJSH1 »BYSM1,IFRR) 

HJPM=-8JMl+DN*ajM/Z 

BJPSM=-SIG*BJSMl+ON*0JSM/Z 

8YPM=-BYMl+0N*BYM/Z 

BYPSM=-S IG4BYSM1+0N+BYSM/Z 

BJ PP'^ = BJ « I /Z- I ON / ( Z* Z } - ( ON ZZ ) **2+ 1 . 00 1 J M 
BYPPM=RYMl/Z-(ON/{ Z*Z )- ( ON /Z ) ♦♦2+1 . 00) *BYM 
BJPPSM=BJSM1^S IG/Z-ION/{ Z^ Z ) - I ON/7 ) ♦♦2+S I G^*2 ) *BJSM 
BYPPSM = BYSMl + S lG/7-( ON / ( Z )- ( DN/Z ) ♦♦? + SI &♦♦? ) ♦BYSM 
BP=flJPPH^HYPSH+B JPM*BYPPSM-R JPPSM^BYPM-BJPSM^BYPPH 
GO TO 81 

80 CALL OBESJYIZ, l,BJi,BYl*BJ2fBY2fIERR) 

ZS=SIG*Z 

CALL OBESJYIZS » 1, R JS 1, BY SI ,B J S2 f BY S2 1 1 ERR ) 

B=SlG^l BJ1*BYS1-BJS1^BY1 ) 

4 
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BP~SIG*BJl*(-S IG*BYS2<^BYSl/7 ) «-S I G^BYSl * (-8 J2«-B J 1/ Z J-SIG* 

ASJ S1*I -RY2 + BY1/Z J -SI G*8 Yl* ( -SIG*B JS2+B JSl/ Z ) 

WRITE{6,4| B,BP 
R1 CONTINUE 

IFtBPJ 61,61,63 
61 INDN=-1 
GO TO 65 
63 INDP=1 
65 IND=INDN«-INOP 

IF { INO.FQ.OJ GO TO 70 
GO TO 60 

17 Z = 3.000 / I l.DO - SIGI 
1=1 

70 CIJ + 1)=Z+. lOO*DFLOAT(n 
30 CONTINUE 
WR ITEI 6, 00) 

90 FORMAT! » 1* ,10X, • ANNULAR DUCT E IGENVALUES* ,//» 

WRITE!6,91I SIG 

91 FORMATUH ,10X,*RA0IUS RATIO =*,020.14//) 

WRITE! 6, 92) HM 

92 FORMATUH , lOX ,' ANGULAR WAVE NUMBER =»,I5) 

00 94 I = 1,NNZ 

NI=I-1 

94 WR ITE! 6,93) NI,ROnT(I) 

93=!F0RMAT( IH ,10X,*RADIAL WAVE NUMBER =• , I 5 * I OX, • F IGE NVALUE 
A KMN =*,025.16//) 

RETURN 

- CALCULATE VELOCITIES - _ _ - 

ENTRY VFLCAL ( lA, BESM) 

NW=IA-l 
XKMN=ROQT( lA) 

IF{ ( MW.EO.O) .AND.! NW.EO. 01 ) GO TO 200 
CALL OBESJYI X< MN , MW, B J ,B Y, BJ J, BYY, lERR) 

8J PM=- BJ J+DFLO AT I MW I *B J / XK HN 
BYPM=-BYV+DFL3AT(MW»*BY/XKMN 

CALCULATE THF MODAL EIGENFUNCTION CONSTANT QMN 
QMN=-BJPM/8YPM 
GO TO 210 
200 QMN=0.000 
210 CONTINUE 

WRITE! 6, 5) QMN 

5 FORMATUH ,10X,*QMN =*,015.8/1 

C CALCULATE MODAL VELOCITY PROFILE ON THE DUCT FACE 

00 50 I = 1,NARG 
IFU-? ) 51,51,52 

51 R! I»=RARGI II*BB/2.D0+ AINK 
GO TO 57 

52 ie=i-i 
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n=R( IE)<-CRARGUei +RARGI I n*BB/2-D0 
C - - THIS CANI^nT BE USED WITH CIRCULAR DUCT DATA - - 

ST IF tINK.EQ.O .AND, R(Il.LT.AA) GO TO 170 

6ARG=R( D+XKMN/BB 

IF(IMW,EQ.O).AND.(NW.EQ.On GD TO 150 
CALL DBESJYI BARGtHW, BESJ , BESY, BUM, RUMM, lERR) 

GO TO 160 
150 BFSJ=1.D0 
BESY^O.DO 
160 CONTINUF 

BESH( I )=BESJ+3MN*BESY 

GO TO 50 

170 BESMI I ) = 0.00 
50 CONTINUE 

WRiTE(6|55) NW 

55 FORMATIIH t lOX , • VELOC I TY DISTRIBUTION ACROSS THE DUCT FACE FOR RAD 
UAL MODE NO = S 15) 

DO 54 J=1,NARG 

54 WRITE(6»53) RARGU ) » R IJ I ,BESH( J I 

53 FORMATIIH ,*R DELTA • , 2X , D1 5. 8, 2X, " R ARGUMENT* ,2X, D15 ,8 »2X V FLOCI T 
AY DISTRIBUTION*, 2X, 015. 8/) 

11 CCNTINUF 
RETURN 
END 


SUBROUTINE DBESJYI X,NO, BJ,BY, BJJ,8YY,IERR> 

C ZDH915 88 - 25 

DOUBLE PRECISION DFLOAT,DLOG 
DOUBLE PRECISION J ( 350) , JO, Jl, JOO 

DOUBLE PRECISION YI 80) , YO, Y1 , X,C,B J,B Y,Pl , RN, 8JJ, BYY 
PI =3. 1 41 592653 589793D0 
lE(X.GT.O.DO) GO TO 90 
IERR=0 
GO TO 80 
90 CONTINUE 
IERR=1 
FAC = 1.5 

IE (X ,GT. 50,0) FAC = .5 

NX=3.9+FAC*X 
JI 2*NX+2) = O,D0 
J( 2*NX<-l)=l«0-29 
C=O.DO 
JO=O.DO 
Jl=0.00 
Y0=0.00 
Y1=0,00 
DO 20 NB=1,NX 
N=NX+1-N8 
RN = DFLOATlN ) 

JI 2*N)=2.D0*I 2,DO*RNl-l,DO)*J(2*N^-l) /X-J(2*N4-2 ) 

J( 2*N-l)=4.D0*RN*J(2*NI/X-JI2*N-H) 

Y0=Y0+I-1)**N*J( 2*N)/PN 
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Yl=YL + {-l)**N*(2*D0*RN*UnO>*J(2*N+ll/RN/(1.0O+RN> 
20 C=C+2,D0*J{ 2*N ) 

JOO=2.DO*J( H/X“J(2I 
C=J00+C 

YO=(DLOG(X/2.Dn I +. 57721 5 66 490153001 *J00-2.O0*Y0 

Yl=fOLOG{X/2.00l-.42?78433509847D0I^Jl ll-JOO/ X-Yl 

YO = 2.OO*Y0/C/P I 

Yl = 2.00*Yl/C/P i 

JO=JOO/C 

V{ U=Yl 

Y( 2»=2.D0*Y( 1) /X-YO 

NNX=2+^G 

DO 30 I = ltNNX 

Y( I+2)=2.D0*0FLnATl !♦!)*¥( I+U/X-YU ) 

30 ji n=jn )/c 

lFtNO.EQ,OI GO TO 70 
BJ =J(MD) 

BY=Y(NO) 

NP=NO+l 
8JJ=JtNP ) 

BYY=Y( NP ) 

GO TO 80 
70 BJ=JO 
BY=YO 
BJJ=J( 1) 

BYY=Yl n 
00 CONTIMUE 
RETURN 
END 
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APPENDIX 8 - EXAMPLE PROBLEMS 


LISTING OF INPUT DATA 


SCMl- INFINITE CIRCULAR DUCT 


^■READ 

PUNCH=T, 

DLU^F, 

+ END 




iOlbAO 

? B 0 0 0 

0 





n.o 

1.0 

343, CO 

343. n 

c, rool 

0,03 


0.02S0 

o.n?5 

0.050 

0.05 

C, 05 0 

0,050 

0.075 

0,n7‘> 

0.075 

0.075 

0. 075 

0. 075 

0,10 

0, 10 

0. 10 

P.IO 

0. 10 

0. 1 

0. 1 

0.1 

0.1 

0. IZb 

0. 125 

0, 125 

0. 125 

0. 125 

0, 125 

0, 125 


0.125 

0.125 

C.1Z5 

0. 125 

0.125 

0.125 

f). 15 

0.15 

0. 150 

0. 150 

Q. 150 



0.125 

0.125 

0 , loco 

0.0750 

,0500000 

0.0250 

0.025 

0.0500 

0.0500 

0. 0750 

n. 0750 

0. 0750 

0.0500 

0.0500 

0.025 

0,025 






01 







l.noon 

l.OOOO 

1.0000 

1. 0000 

1. OOOO 

l.onon 

1,0000 

1.0000 

1 . 0000 

l.OCCO 

1,0000 

1, OOOG 

1,0000 

1.0000 

l.oooo 

l.oooo 







ENC 







SFrtl- 

INFINITF ANNULAR DUCT 





+ RFAI) 

PUNCH=T , 

GL ()=F, 

+ END 




301540 

2 B n 0 0 

0 





0.0 

1.0 

3/^3. 00 

343, 0 

C. 0001 

0.03 


0.0250 

0,025 

0,050 

C.C5 

C. C50 

0.050 

0.0/5 

0,075 

0.075 

0,075 

0,075 

0. 075 

0,10 

o. 10 

0, 10 

0. 10 

0. 10 

0. 1 

0. 1 

0.1 

0,1 

0.125 

0. 125 

0. 125 

C, 125 

0. 125 

0.125 

0,125 

0.125 

0.125 

0.125 

0. 125 

0, 125 

0. 125 

0.125 

0.15 

0. 15 

0. 150 

0. 150 

n, 150 



0.125 

0. 125 

0 . loco 

0, 0750 

,0500000 

0.0250 

0.025 

0.0500 

0,0500 

0.0750 

G, 075^ 

0, 0750 

0,050n 

C.05O0 

0.025 

0.025 






01 







o.ooon 

0.0000 

0.0000 

G.OOOO 

0. 0100 

0.0000 

1,0000 

1 .Of^OO 

l.oooo 

l.OOCO 

l.oooo 

1. OOOO 

1.0000 

1.0000 

l.nc-oo 

l.OOCO 







ENH 
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SFM !- 

INF, ANN. DUCT VflTH 

CENTFRBODY 




+PF&H 

301010 

PUNCH=T, 

? B 0 1 no 

m.D=F, 


♦ END 



n,*i 

1.0 

343.no 

343. 0 

0.0001 

0.03 


• 0?F 

• 2 

• 05 
.225 

• O?*: 
.25 

• 1 

.125 

.15 

• 175 

.025 

. 2 

.05 
. 225 

• 075 
.25 

• 1 

.125 

.15 

.175 

.025 

.05 

01 

.05 

.025 

,05 

.025 

.075 

.075 

.075 

.05 

1 .0000 
1.0000 

l.OOOO 

l.OOOC 

1. COCO 

1 . cnoo 

1, 0000 

t.oonr 

1,0000 

1 .noon 


FND 


ANN. DUCT v*ITH CENTFRftODY 
4-prM) PIJMCM^T, 0U)=F, +FNO 

TOK'^0 8 0 1 0?0 


n.b 

1.0 

343.00 

343.0 

0. COOL 

0.03 


.025 

.n^5 

.05 

.0 75 

,075 

.075 

.075 

,n 75 

. 1 

.125 

.15 

.1 75 

.175 

.2 

.2 

.2 

.? 

.2 

.2 

.2 


.025 

,025 

.05 

.075 

.075 

.075 

.075 

.075 

. 1 

.125 

.15 

.1 75 

. 175 

.2 

.2 

. 2 

.2 

.2 

.2 

.2 


.025 

.05 

01 

.025 

.025 

.05 

.025 

.075 

.075 

,075 

.075 

1 .OOOO 

1 .Ofi oo 

1.0000 

1 .0000 

1.0000 

l.OrCO 

l.OOOO 

l.OOOO 

l.OOOO 

1 .or no 


END 
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APPENDiX 8 - EXAMPLE PROBLEMS 
OUTPUT FOR EXAMPLE PROBLEM 1 


PRJBCEDING PAGE BLANK NOT FILMEO 
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BOX METHOD INTEGRATION FOR SOURCE DISTRIBUTION OVER DUCT 


SEMI-INFINITE CIRCULAR DUCT 


DUCT OR = l.OOCO FT DUCT IR = C.O 

FREQUENCY = 343,0000 RAD/SEC SOUND SPEED = 

MODE NUMBER = 0 


FT 

343,0000 


FT/SFC 


+PPNT1 
L = 

+ FND 
+ PRNT? 
N? = 

+ END 


30, M = 


16, N = 


4C,NN= 


2 ,LP = 


0, INK = 


0,CHK= 0,9g<99g9<;0E-0A,XANG= 0, 29999997F-01 


BLOCK SIZE IN THE -Z DIRECTION NORMALIZED WITH RESPECT TO B. 


0.25C0E-01 
0,50C0E-01 
0. 75CQF-01 


0, ICCOE 
0, lOCOF 
0. ICCOE 
0, 1250E 
0. 125CF 
0, 1260F 
0, I 5CGE 


0,25 
0, 50 
0, 75 
0.10 
0 . 10 
0 . 1 ? 
0.12 
0 , 1 ? 
0.12 
0.15 


OOF-01 
OOE-01 
OOE-01 
noE 00 


OOF 00 
50E 00 


50F 00 
50E 00 


50F 00 
OOF 00 


0. 50C0E-01 
0. 7500E-01 
0. 75COE-01 
OclOOOF 00 
0, lOCOE 00 
0, 1250F OC 
0, I250E 00 
0.1250E 00 
0. 1250F 00 
0. 15C0E 30 


0. 5000 
0,7500 
0, 7500 
0, 1000 
0. 1000 
0. 1?50 
0. 1 2 50’ 
0. 1250 
0. 1 500 
C. 1500 


F-ni 
F-OI 
F-f^l 
F 00 


F 00 
F no 


F 00 
F 00 


F 00 
E 00 


BLOCK 


SIZF IN THE 
0..12 50F 00 
0, 5CC0E-01 
0, 5CC0E-01 
0. 5CC0E-01 


DIRECTION NORMALIZED WITH RESPFCT TO B. 

0. 12 50E 00 0. lOOOE 00 0, 7500E-01 

0.25C0E-01 C, 2500F-QI 0, 50n0F-01 

0.7500F-01 0.75C0E-01 0,7500F-01 

0.5000F-01 0. 2500E-01 0, 2500E-01 


BLOCK NO 1 


RADIUS = 


C.062 5 


1 

6.2332 

+ 

O.f 

I 

2 

0.0 

+ 

0,0 

I 

3 

C.O 

+ 

o.c 

I 

4 

0.0 

+ 

o.c 

I 

5 

0.0 

4 - 

0,0 

I 

6 

0.0 

+ 

o.c 

I 

7 

C.O 

+ 

0,0 

I 

3 

0,0 

+ 

o.c 

I 

9 

0.0 

+ 

o.c 

I 

10 

C.O 

+ 

o.c 

I 

11 

0,0 

+ 

o.c 

I 

12 

0,0 

+ 

O.G 

I 

13 

C.O 

+ 

o.c 

I 

14 

C.O 

+ 

o.c 

I 

15 

0.0 

4 - 

o.c 

I 










0 

0.0 

+ 

O.C 

I 



17 

0.0027 


0,0006 

I 

Z = 

-0.0125 

n 

0.0082 

+ 

0.0018 

I 

7 . = 

-0.0375 


0.0325 

+ 

0.C071 

I 

1 = 

-0. 0750 

^ 20 

0.0535 

+ 

0.0118 

I 

1 = 

-0. 1250 


0.0734 

+ 

0.G165 

I 

Z = 

-0. 1750 

n 22 

0.0920 

*■ 

0.0212 

I 

1 = 

-0. 2750 

Ll 23 

0.1694 

+ 

0.04Q4 

I 

z = 

-0.2875 

24 

0.2016 

+ 

0.0507 

I 

z = 

-G. 362 5 

n 25 

0.2273 

♦ 

0.0608 

I 

z = 

-0.4375 

U 26 

0.2466 

+ 

0.0707 

I 

z = 

-0.5125 

27 

0.2600 


0.0804 

I 

1 = 

-0, 5875 

f 1 28 

0.2681 

+ 

0,0898 

I 

z = 

-0. 6625 

29 

0.3624 

+ 

0.1338 

T 

z = 

-0.7500 

^ 30 

0,3601 

4- 

0. 1491 

I 

z = 

-0. 8500 

31 

0.3512 

+ 

0.1635 

I 

z = 

-0. 9500 

D 52 

0, 3374 

4- 

0. 1 769 

I 

7 = 

-1,0500 

U 33 

0.320? 

+ 

0.1892 

I 

/ = 

-1, 1500 

34 

C.3004 

4- 

0.2004 

I 

Z = 

-1.2500 

n 

0.2790 

4 

0.2104 

I 

Z = 

-1. 3500 

U 36 

0.2565 

4 - 

0.2192 

I 

1 = 

-1,4500 


0.2334 

4 - 

0.2266 

I 

/ = 

-1. 5500 

f. 38 

0.2587 

4- 

0.2917 

I 

z = 

-1.6625 

1 [| 39 

0*2219 

4- 

0.2986 

I 

z = 

-1.7875 

L] 40 

0. 1854 

4 

0,3029 

I 

z = 

-lo 9125 

^1 

0.1496 

4- 

0.3044 

I 

z = 

-2. 0375 

1 n '‘2 

0.1147 

4- 

0.3033 

I 

z = 

-2. 1625 


0,0309 

4- 

0,2995 

I 

z = 

-2.2875 

1 ^^4 

0.0486 

4- 

0.2932 

I 

z = 

-2.4125 

n 

0,0178 

4- 

0,2846 

I 

z = 

-2, 5375 


-0.0112 

4- 

0.2736 

I 

7 = 

-2,6625 

i ^ 47 

-0,0382 

+ 

0,2605 

I 

1 = 

-2.7875 

1 ... 

-0.0632 

♦ 

0.2454 

I 

7 = 

-2. 9125 

49 

-0.0859 

4- 

0,2286 

I 

Z = 

0375 

U 50 

-C. 1063 

+ 

0.2103 

I 

Z = 

-3. 1625 

1 31 

-0.1243 

+ 

0,1906 

I 

1 = 

-3.2875 

n 32 

-0,1694 

4^ 

0,2013 

I 

7 = 

-3,4250 

U 33 

-0.1873 

4- 

0,1 700 

I 

Z = 

-3. 5750 

54 

-0.2008 

4- 

0.1379 

I 

2 = 

-3. 7250 

n 35 

-0.2098 

4- 

0.1054 

I 

Z = 

-3.8750 

0 

-0,2145 

4- 

0.0731 

I 

Z = 

-4.0250 

jj BLOCK NO 

2 RADIUS = 

r.1875 



n 1 

0.0 

♦ 

0,0 

I 



U 2 

6.2832 

4- 

O.C 

I 



Li .j 

O.C 

4- 

O.C 

I 



4 

C.O 

4- 

O.C 

I 



0 3 

0,0 ■ 

4- 

O.C 

1 



0 6 

0.0 

4^ 

O.C 

I 



f 7 

0,0 

4- 

0,0 

I 



I 1 
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1 ^ 








8 

0.0 

•f 

o.c 

I 



9 

0.0 

♦ 

0.0 

1 



10 

0*0 

+ 

0.0 

I 



11 

0.0 

* 

0.0 

I 



12 

0.0 

* 

0.0 

I 



13 

0.0 

* 

0.0 

I 



1^' 

0.0 


0.0 

I 



I'i 

0.0 

+ 

0.0 

I 



16 

0.0 

+ 

0.0 

I 



17 

0.0029 

+ 

0.0006 

I 

z = 

-0.0X25 

IB 

0.0086 

+ 

0.0018 

I 

z = 

-0.0375 

19 

0.0343 

+ 

0.0071 

I 

z = 

-0.0750 

20 

0.0563 

+ 

0.01,18 

I 

z = 

-0, 1250 

21 

0.077? 

+ 

0.C164 

I 

z = 

-0. 1750 

22 

0^0965 

* 

0.0211 

I 

7 = 

-0.2250 

23 

0. 1768 

+ . 

0.C4O3 

I 

z = 

-0.2875 

24 

0.2092 


0,0506 

I 

z = 

-n, 362 5 

25 

0.2345 

+ 

0.C607 

I 

z = 

-0.4375 

26 

0..2530 

* 

0,0705 

I 

2 = 

-0.5125 

27 

0.2653 

+ 

0,0802 

I 

/ = 

-0.5875 

28 

0.2723 

+ 

0.0895 

I 

z = 

-0.6625 

29 

0.3664 

+ 

0.1354 

I 

/ = 

-0.7500 

30 

0.3624 

+ 

0.1486 

I 

z = 

-0. 8500 

31 

0.3522 

+ 

0.1630 

I 

z = 

-0. 9500 

32 

0.3375 

♦ 

0,1764 

1 

2 = 

-1.0500 

33 

0-3195 

♦ 

0.1887 

I 

z = 

-1. 1500 

34 

0.2993 

+ 

0.1998 

I 

z = 

-1.2500 

35 

0.2777 

- 1 - 

0.2098 

I 

z = 

-1.3500 

36 

C.2550 

*■ 

0.2185 

I 

z = 

-1.4500 

37 

C. 2318 


0.2259 

I 

z = 

-1.5500 

38 

0.2560 

+ 

0.2908 

I 

1 - 

-1,6625 

39 

0.2200 

+ 

0.2977 

I 

1 = 

-1.7875 

40 

0.1837 

+ 

0.3018 

1 

z = 

-1.9125 

41 

0. 1480 

♦ 

0.3033 

T 

z = 

-2.0375 

42 

0.1132 

+ 

0.3022 

I 

z = 

-2. 1625 

43 

0.C796 


0.2984 

I 

z = 

-2.2875 

44 

0.0474 

+ 

0,2922 

I 

z = 

-2,412 5 

45 

0.0168 

+ 

0.2835 

I 

z = 

-2,5375 

46 

-0.0120 

+ 

0.2725 

I 

z = 

-2.6625 

47 

-0.0390 


0.2594 

I 

z = 

-2.7875 

48 

-0,0638 

+ 

0.2444 

I 

z = 

-2,9125 

49 

-0.0864 

♦ 

0,2276 

I 

2 = 

“3,0375 

50 

-0. 1067 

+ 

0.2093 

I 

z = 

-3. 1625 

51 

-0.1246 

+ 

0,1897 

I 

z = 

-3.2875 

52 

-0. 1696 


0.2002 

I 

z = 

-3,4250 

53 

-0.1874 

+ 

0,1690 

I 

z = 

-3.5750 

54 

-0.2008 

+ 

0,1369 

I 

1 = 

-3.7250 

55 

-0. 2097 

+ 

0.1045 

I 

z = 

-3. 8750 

56 

-C.2143 

♦ 

0.0723 

I 

1 = 

-4.0250 
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BLOCK NO 3 RADIUS = C,3000 


1 

0,0 

•I- 

0.0 

I 



2 

0.0 

+ 

o.c 

I 



1 

6.2832 

+ 

0.0 

I 



4 

0.0 

*■ 

0,0 

I 



5 

0.0 


o.c 

I 



6 

0.0 

+ 

o.c 

I 



7 

0.0 

+ 

0.0 

I 



R 

0.0 


0.0 

I 



q 

0.0 

+ 

0.0 

I 



10 

0.0 

+ 

0,0 

I 



u 

c.o 

+ 

0.0 

I 



12 

0.0 

+ 

0.0 

I 



13 

0.0 

+ 

0,0 

I 



14 

0.0 

+ 

0,0 

I 



15 

0.0 

* 

o.c 

I 



16 

c.o 

+ 

0.0 

I 



17 

0.0032 

+ 

0.0006 

I 

2 = 

-0.0125 

18 

0.0096 


O.C018 

I 

7 = 

-0, 0375 

IR 

0.0379 

+ 

0.007C 

I 

Z = 

-n. 0750 

20 

0.0620 

+ 

o.c 11 7 

I 

Z = 

-0. 1250 

21 

0.0846 

-1- 

0.C164 

1 

Z = 

-0, 1750 

22 

0. 1051 

+ 

0.C210 

I 

Z = 

-C.2250 

23 

0.1910 

+ 

0.C401 

I 

7 = 

-0.2875 

24 

C.2235 

+ 

0.C503 

I 

Z = 

-0, 362 5 

25 

0.2478 

<1- 

O.C 60 3 

I 

7 = 

-0.4375 

26 

0.2645 

+ 

0.0702 

I 

Z ~ 

-0. 5125 

27 

0. 2747 

+ 

0,0797 

I 

7 = 

--0. 5R75 

28 

0.2796 

+ 

0.C890 

I 

Z - 

“0.662 5 

29 

0.3729 

+ 

0.1327 

I 

Z = 

-0,7500 

30 

0,3661 

+ 

0.1478 

I 

Z = 

-0.8500 

31 

0,3536 

+ 

0.1621 

I 

z = 

-0,9500 

32 

0.3373 

4 

0.1754 

I 

z = 

-1,0500 

33 

0.3182 

4 - 

0,1876 

I 

z = 

-1. 1500 

34 

0,2973 

4 - 

0.1987 

I 

z = 

-1,2500 

35 

0.2752 

+ 

0,2086 

I 

z = 

-1.3500' 

36 

0.2523 

4 - 

0,2173 

I 

z = 

-1.4500 

37 

0,2290 

4 - 

0,2246 

I 

7 = 

-1.5500 

38 

0.2533 

4 

0,2891 

I 

z = 

-1.6625 

39 

0.2167 

4 - 

0,2959 

I 

z = 

-1, 7875 

40 

0.1806 

4 - 

0.3001 

I 

2 = 

-1.9125 

41 

0,1451 

4 - 

0.3015 

I 

Z = 

-2.0375 

42 

0,1106 

4 - 

0,3003 

I 

7 = 

-2.1625 

43 

0.0773 

* 

0.2966 

I 

1 = 

-2.2875 

44 

0.0454 

4 

0.2903 

I 

7 = 

-2.4125 

45 

0.C15C 

4- 

0.2816 

I 

2 = 

-2.5375 

46 

-0.0136 

4 - 

0,2706 

I 

Z = 

-2.6625 

47 

-0.0402 

4- 

0.2576 

I 

z = 

-2. 7875 

48 

-0,0649 

+ 

0.2426 

I 

z = 

-2.9125 

49 

-0.0873 

4 - 

0,2258 

I 

z = 

-3.0375 

50 

-0.1074 

+ 

0.2076 

I 

z = 

-3,1625 


51 

-0.1251 

♦ 

0.1880 

I 

Z = 

-3.2875 

52 

“0. 1700 

♦ 

0.1983 

1 

z = 

-3.4250 

53 

-0. 1875 

+ 

0,1672 

I 

z = 

-3.5750 

5^ 

-0.2007 

+ 

0.1353 

I 

z = 

-3.7250 

55 

-0.2095 

+ 

0.1030 

I 

1 = 

-3.8750 

56 

-0.2139 

+ 

0.0709 

I 

z = 

-4,0250 

BLOCK NO 

A RADIUS = 

0.3875 





1 

0.0 

+ 

0.0 

I 



2 

0.0 

+ 

0.0 

I 



3 

0.0 

4 

O.G 

I 



4 

6.2832 

4 

O.G 

1 



5 

0.0 

4 

0,0 

I 



6 

0.0 

4 

0,0 

I 



7 

0.0 

4 

0.0 

I 



8 

0.0 

4 

0,0 

I 



9 

0.0 

4 

0,0 

I 



10 

c.o 

4 

0.0 

I 



11 

0.0 

4 

0,0 

I 



12 

0.0 

4 

0.0 

I 



13 

0.0 

4 

o.c 

I 



14 
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Z = 

-1. 0500 

33 

-0.0022 


O.G818 

I 

z = 

-1. 1500 

34 

0. 0040 

♦ 

0. C868 

I 

z = 

-1. 2500 

35 

0.0109 

+ 

0.0919 

I 

7 = 

-1.3500 

36 

0.0184 

+ 

0.0968 

I 

Z = 

-1,4500 

37 

0. 0266 

+ 

0.1017 

I 

Z = 

-1, 5500 

33 

0. C460 

+ 

0,1338 

I 

Z = 

-1. 6625 

39 

0.0616 


0.1409 

I 

z - 

-1. 7875 

40 

0. 0788 

+ 

0.1478 

I 

Z = 

-1.9125 

41 

0,0978 

+ 

0.1 543 

I 

2 = 

-2.0375 

42 

0. 1187 

+ 

0^1 604 

I 

1 - 

-2, 1625 

43 

0. 1417 

+ 

0.1660 

T 

1 = 

-2, 2875 

44 

0- 1669 

+ 

0.1 712 

I 

7 = 

-2.4125 

4b 

0. 1946 

+ 

0.1759 

I 

1 - 

-2. 5375 

46 

0.2252 

+ 

0,1800 

I 

1 = 

-2. 6625 

47 

0. 2593 

4^ 

0,1836 

I 

Z - 

-2. 7875 

48 

0.2981 

+ 

0.1865 

I 

z = 

-2. 9 1 5 

49 

0.3436 

+ 

0.1888 

I 

z = 

-3. 0375 

50 

0.401C 

+ 

0.1905 

I 

1 = 

-3. 1625 

51 

0.4820 

4 

0,1916 

I 

( = 

-3, 2875 

52 

C.8630 

4 

0.2304 

I 

7 = 

-3,4250 

53 

0.5676 

4 

0.2298 

I 

Z = 

“3. 5 75 0 

54 

0.4588 

4 

0.2281 

I 

1 = 

-3, 7250 

55 

0, 3839 

4- 

0,2253 

I 

Z = 

-3.8 75 0 

56 

0.3246 

4 

0,2214 

I 

z = 

-4,0250 

BLOCK 

NO 53 Z DISTANCE 

= -3. 

5750 




1 

-0.0027 ♦ 

O.C025 I 

P = 

0.0625 

2 

-0.0082 + 

0.0073 I 

R = 

0. 187 5 
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3 

-C.0105 

+ 

0.C092 

I 

P = 

0. 3000 

4 

-0#0102 

+ 

0,0087 

I 

P = 

0, 3875 

5 

-0.0079 

+ 

0.C066 

I 

P = 

C.450n 

5 

-C.0043 

+ 

0,0035 

I 

R = 

0. 487 5 

7 

-0.0045 

+ 

0,0036 

I 

P = 

n.5l2 5 


-0. 0097 

+ 

0,C077 

1 

R = 

0.5500' 

9 

-C.0106 

+ 

0.C082 

I 

P = 

0,6000 

10 

-0.0177 

•i- 

0.C130 

I 

R = 

0. 6625 

11 

-0.0197 

+ 

0.0138 

I 

P = 

0.7375 

12 

-0.0218 

+ 

0,014A 

1 

R = 

0, 8125 

13 

-C.0157 

f 

0.C098 

I 

R = 

0.8750 

14 

-0.0166 

+ 

0.0099 

I 

R = 

0,9250 

15 

-C.0086 

4 

0.C049 

I 

P = 

C.9625 

16 

-0.0089 

4 

0.0049 

I 

P = 

0.9875 

17 

-0.0090 

4 

0.CQ51 

I 

Z = 

-0.0125 

18 

-0,0090 

4 

0,0053 

I 

L = 

-0. 0375 

19 

-C.0179 

4 

0.0114 

I 

2 = 

-0.0750 

20 

-0.0178 

4 

0,0124 

I 

Z == 

-0, 1250 

21 

-C.0177 

4 

0.C135 

I 

7 = 

-0. I 750 

2? 

“0.0175 

4 

0.0146 

I 

1 = 

-0. 2250 

23 

-0.0259 

4 

0.C239 

I 

Z = 

-C. 2875 

24 

-C.0253 

4 

0.C265 

I 

7 

-0, 3625 

25 

-0.0245 

4 

0.0290 

I 

z = 

-0,4375 

26 

-C.0236 

4 

0.C31 7 

I 

7 = 

-0, 5125 

27 

-0.0224 

4 

0.C344 

I 

Z = 

-0, 5875 

28 

-0,021 1 

4 

0.0371 

I 

7 = 

-0. 662 5 

29 

-C, 0258 

4 

O.C53tt 

I 

7 = 

-0, 7500 

30 

-0.0227 

4 

O.C588 

I 

7 = 

-0. 8500 

31 

-0.0191 

4 

0.C639 

I 

Z = 

-C.9500 

32 

-0.0150 

4 

0.0690 

I 

Z = 

-1.0500 

33 

-0.0103 

4 

0.0741 

I 

7 = 

-1. 1500 

34 

-0,0050 

4 

0.C792 

I 

7 = 

-1. 2500 

35 

G.C009 

4 

0.C843 

I 

7 = 

-1. 3500 

36 

0.0074 

4 

0,0894 

I 

7 = 

-1.4500 

37 

0.0145 

4 

0.C944 

I 

Z = 

-1. 5500 

38 

0.0293 

4 

0.1248 

I 

z = 

-1. 6625 

39 

0.0430 

4 

0,1323 

I 

Z = 

-1. 7875 

40 

0.0583 

4 

0.1395 

I 

7 = 

-1.9125 

41 

0.0752 

4 

0,1464 

T 

Z = 

-2. 0375 

42 

0.0939 

4 

0. 1 530 

I 

z = 

-2, 1625 

43 

0.1144 

4 

0.1592 

I 

7 = 

-2, 2875 

44 

0. 1369 

4 

0.1649 

T 

Z = 

“2,4125 

45 

0. 1617 

4 

0.1702 

I 

1 = 

-2. 5375 

46 

0,1989 

4 

0,1750 

I 

Z = 

-2.6625 

4 7 

0. 21 89 

4 

0.1792 

I 

7 = 

-2. 7875 

4H 

0.2522 

4 

0.1829 

I 

Z = 

-2.Q125 

49 

0.2899 

4 

0,1860 

I 

1 = 

-3.0375 

50 

C. 3338 

4 

0. 1 884 

I 

1 = 

-3. 1625 

51 

C.3881 

4 

0.1903 

I 

Z = 

-3,2875 

52 

0. 5676 

4 

0.2298 

I 

7 = 

-3,4250 

53 

0, 8630 

4 

0,2304 

I 

Z = 

-3, 5 75 0 

54 

0, 5676 

4 

0.2298 

T 

Z = 

-3.7250 

55 

0.4588 

4 

0.2201 

I 

Z = 

-3. 8750 
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0.3839 + 0.2283 [ 

Z = 

-4. 0250 

BLOCK 

NO 5A 2 niSTANCF- = -3.7250 




1 

-0.0028 

♦ 

0.0019 

I 

R = 

0.0625 

2 

-0.0084 

+ 

0.0057 

I 

H = 

0. 1875 

3 

-0.0108 

+ 

0.C071 

I 

P = 

'0, 3000 

4 

-0.0105 

+ 

0.C067 

I 

P = 

0.38T5 

5 

-0.0081 

♦ 

0.0050 

I 

R = 

0.45P0 

6 

-0.0044 


0.0027 

I 

R = 

0« 487 5 

7 

-0.0046 

♦ 

0.0028 

I 

P = 

0. 512 5 

B 

-C.0099 

+ 

O.OC59 

I 

R = 

0, 5500 

9 

-0.0108 

+ 

0.0062 

I 

P = 

0. 600 0 

10 

-0.0179 

♦ 

0.0098 

T 

R = 

0. 6625 

11 

-0.0199 

+ 

0.0103 

I 

P = 

0.7375 

12 

-0.0218 

-f 

0,0106 

I 

P = 

0. 812 5 

13 

-0.0156 

+ 

0.C071 

I 

R = 

0. 8 75 0 

14 

-0,0165 

+ 

0.C071 

I 

P = 

0.9250 

15 

-0.0086 

+ 

0.C035 

I 

P = 

0. 962 5 

16 

-0,0088 

+ 

0.C035 

I 

R = 

0.9875 

17 

-0.0089 

+ 

0.CC36 

I 

7 = 

-0.012 5 

18 

-0,0089 

+ 

0.0038 

I 

7 = 

-0. 0375 

19 

-C.0179 

+ 

0.C084 

I 

7 = 

-0. 0750 

20 

-0.0179 

+ 

0.0094 

I 

7 = 

-0. 12 50 

21 

-0.0179 


0,0104 

I 

Z = 

-0. 1750 

22 

-0,0179 

+ 

0.0114 

I 

7 = 

-0. 2250 

23 

-C.0267 

+ 

0.0190 

I 

Z = 

-0, 2 875 

24 

-0, 0263 

♦ 

0,0214 

I 

7 = 

-0, 362 5 

25 

-C.0259 

+ 

0,0239 

I 

7 = 

-0. 4375 

26 

-0,0253 


0.0265 

I 

Z = 

-0,5125 

27 

-0.0245 

+ 

0,0290 

I 

7 = 

-0. 5875 

28 

-0.0236 


0.C31 7 

I 

Z = 

-0.6625 

29 

-C. 0296 

+ 

0.C465 

I 

Z = 

-0. 7500 

30 

-0,0272 

♦ 

0.0513 

I 

z = 

-0, 8500 

31 

-0,0243 

+ 

0.0563 

1 

7 = 

-0.9500 

32 

-0.0210 

+ 

0.C613 

I 

2 = 

~1, 0500 

33 

-0,0171 

+ 

0.C664 

I 

7 = 

-1, 1500 

34 

-0.0127 

+ 

0.C71 5 

r 

7 = 

-1,2500 

35 

-0.0077 

+ 

0.C767 

I 

7 = 

-1.3500 

36 

-0.Q022 

+ 

G.cei e 

I 

7 = 

-1.4500 

37 

0,0040 

* 

0.C868 

I 

2 = 

-1, 5500 

39 

0.0147 

♦ 

0.1156 

I 

7 = 

-1.6625 

39 

0,0267 


0.1233 

I 

7 = 

-1.7875 

40 

0,0402 

+ 

0.1308 

I 

7 = 

-1.9125 

41 

0,0551 

+ 

0.1381 

I 

2 = 

-2, 0375 

42 

0.0717 

+ 

0.1451 

I 

Z = 

-2. 1625 

43 

0. 0900 

+ 

0,1517 

i 

2 = 

-2. 2875 

44 

C. 1101 

*■ 

0.1 580 

I 

2 = 

-?.41 2 5 

45 

0.1323 

+ 

0.1638 

I 

2 = 

-2. 5375 

46 

0.1565 

+ 

0.1692 

I 

7 = 

-2. 6625 

47 

C.1832 

+ 

0,1741 

I 

2 = 

-2. 7875 

48 

0,2126 

+ 

0.1 784 

I 

7 = 

-2. 9125 
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4 





49 

0.2452 

+ 

0. 1822 

I 

1 

=: 

-3,0375 

50 

0,2819 

+ 

0,1854 

I 

1 

= 

-3. 1625 

51 

0,3243 

+ 

0.1880 

I 

Z 

= 

-3. 2875 

52 

0.4588 

+ 

0.2281 

I 

z 


-3.4250 

53 

0.5676 

4 - 

0,2298 

I 

z 

= 

-3. 5750 

54 

0. 863C 

+ 

0.2304 

I 

z 

= 

-3. 7250 

55 

0. 5676 


0,2298 

I 

z 

= 

-3. 8750 

56 

0,4588 

+ 

0.2281 

I 

z 


-4, 0250 

BLOCK 

NO 55 Z 

riSTANCE 

= -3. 

8750 





1 

-0.0028 

+ 

0.0014 

I 

K = 

0. 0625 

2 

-0.0085 

+ 

0.0C42 

I 

P = 

0, 1875 

3 

-0,0108 

+ 

0.0052 

I 

P = 

0. 3000 

4 

-0.01G5 

+ 

0.0048 

I 

R = 

0. 3875 

5 

-c.ooai 

4 

0.0036 

I 

R = 

0.4500 

6 

-0.0044 

4 

0.C019 

I 

P = 

0.4875 

r 

-0.0046 

+ 

0.0020 

I 

P = 

Oo 512 5 

8 

-0.0098 

4 

0.0042 

I 

R = 

0, 5500 

9 

-0.0107 

4 

0.C044 

I 

P = 

0. 6000 

10 

-0.0177 

4 

0.C068 

I 

P = 

0. 6625 

11 

-0.0196 

4 

0.0071 

I 

P = 

0. 7375 

12 

-C.0215 

4 

0.0071 

I 

P = 

0, 8 12 5 

13 

-0, 0154 

4 

0.0047 

I 

P = 

0. 8750 

14 

-0.0162 

4 

0,0046 

I 

P = 

0, 9250 

15 

-0.0084 

4 

0,0022 

I 

P = 

0. 962 5 

16 

—0.0086 

4 

0.0022 

1 

P = 

0. 9875 

17 

-O.OOP7 

4 

0.C023 

T 

z = 

-0.0125 

18 

-0.0087 

4 

0,0025 

I 

z = 

-0. 0375 

19 

-0.0176 

4 

0.0056 

I 

z = 

-0, 075 0 

20 

-0,0177 

4 

O.C065 

I 

z = 

-0, 12 50 

21 

-0.0178 

4 

0.0075 

f 

z = 

-n. 1 750 

22 

-0.0179 

4 

0,0084 

j 

7 - 

-0. 2250 

23 

-0.0269 

4 

0.0144 

I 

Z = 

-0. 2875 

24 

-0.0268 

4 

0,.C167 

I 

1 = 

-0. 362 5 

25 

-C.0267 

4 

0.C19C 

I 

Z = 

-0.4375 

26 

-0.0263 

4 

0.C214 

I 

Z = 

-0. 5125 

27 

-0.025° 

4 

0.0239 

I 

Z = 

-0. 5875 

28 

-0.0253 

4 

0,0265 

I 

Z = 

-0. 662 5 

29 

-0.0325 

4 

0.0393 

I 

7 = 

-0, 7500 

30 

-0.0307 

4 

0.0440 

t 

Z = 

-0. 8500 

31 

-0.0285 

4 

0.0489 

I 

Z = 

-0. 9500 

32 

-C.0258 

4 

0.C538 

I 

7 = 

-1, 0500 

33 

-0.0227 

4 

0.0588 

1 

7 = 

-1. 1500 

34 

-0.0191 

4 

0.C639 

t 

Z = 

-1. 2500 

35 

-0,0150 

4 

0,0690 

I 

Z = 

-1. 3500 

36 

-0.0103 

4 

0.0741 

1 

1 = 

-1. A5C0' 

37 

-0.0050 

4 

0.0792 

I 

2 = 

-1. 5500 

38 

0.0021 

4 

0.1062 

I 

7 = 

-1.662 5' 

39 

0.0125 

4 

0,1141 

I 

7 = 

. -1,7875 

40 

C.0242 

4 

0.1218 

1 

7 = 

-1,9125 


191 


41 

0,0374 

* 

0. 1293 

I 

7 


-2.0375 

42 

0.0520 

+ 

0.1367 

I 

1 


-2, 1625 

43 

0.0683 

*■ 

0.1437 

I 

I 

= 

-2. 2875 

44 

0.0862 

+ 

0,1504 

I 

1 


-2.4125 

45 

0* 1059 


0.1568 

I 

1 

= 

-2, 5375 

4S 

0. 1277 

* 

0,1627 

I 

7 


-2, 6625 

47 

0.1515 

♦ 

0.1682 

I 

1 

= 

-2.7875 

4H 

0.1777 

♦ 

0.1731 

I 

1 

= 

-2.9125 

49 

0.2065 

*■ 

0.1776 

I 

7 

= 

-3,0375 

50 

0.2384 


0,1 81 5 

I 

1 

= 

-3, 1625 

51 

0.2742 

+ 

0.1 848 

I 

1 

= 

-3.2B75 

52 

0,3839 

+ 

0.2253 

I 

1 

- 

-3.4250 

53 

0.4588 

+ 

0,2281 

I 

1 

= 

-3. 5 75 0 

54 

0.5676 

+ 

0.2298 

I 

2 

- 

-3, 7250 

55 

C.8630 

♦ 

0.2304 

1 

1 

= 

-3.8750 

56 

0.5676 

+ 

0.2298 

I 

I 

= 

-4,0250 


Block 

NO 56 Z CISTANCE 

= -4, 

0250 




1 

-0,0028 

-t- 

0.0009 

I 

P 

= 

0.062 5 

2 

-0,0083 

4- 

0.C028 

I 

R 

= 

0. 1875 

3 

-0,0106 

+ 

0.CC34 

I 

R 

= 

0. 300C 

4 

-0,0102 


0.0031 

I 

P 

= 

0.3875 

5 

-0.0079 

+ 

0,C023 

I 

R 

= 

0.4500 

6 

-0.0043 

+ 

0,0012 

I 

R 

= 

0.4875 

7 

-0,0045 

+ 

O.CCl 3 

I 

P 

= 

0.5125 

8 

-0.0096 

+ 

0.0026 

I 

R 

- 

0, 5500 

9 

-0.0104 

♦ 

0.0027 

I 

R 

= 

0. 6nC0 

10 

-0.0172 

+ 

0,C041 

I 

R 

= 

0.6625 

11 

-0^0190 

♦ 

0,0041 

I 

P 

= 

G. 7375 

12 

-0,0208 

♦ 

0.0039 

I 

R 


0. 8125 

13 

-0,0148 

+ 

0.C025 

I 

P 


0. 876 0 

14 

-C.C155 

+ 

0,0023 

I 

P 

= 

0.9250 

15 

-0.0080 

+ 

o.roi 0 

I 

R 

= 

0, 9625 

16 

-0,0082 

♦ 

0,0010 

I 

P 

= 

0.9875 

17 

“0,0083 

« 

0,0010 

I 

Z 

= 

-0.0125 

18 

-C.0084 

+ 

0,0012 

I 

2 

= 

-0.0375 

n 

-0.017C 

♦ 

0,0031 

I 

1 

= 

-0.0750 

20 

-0.0172 


0.0039 

I 

1 

= 

-0. 1250 

21 

-0.01 74 

♦ 

0,C046 

I 

7 


-0. 1750 

22 

-0.0176 

+ 

0,0056 

I 

7 

= 

-0. 2250 

23 

-0.0266 

+ 

O.CIOI 

I 

7 


-0.2875 

24 

-0.0268 

+ 

0.C122 

I 

L 

= 

-r, 362 5 

25 

-0.0269 

+ 

0.0144 

I 

7 

= 

-0,4375 

26 

-0.026P 


0.C167 

I 

1 


-0.5125 

27 

-0.0267 

+ 

0,0190 

T 

7 

- 

-C.5875 

28 

-0.0263 

♦ 

0.0214 

I 

1 


-0. 66? 5 

29 

-C.0344 

+ 

0,0324 

I 

1 


-0. 7500 

30 

-0,0332 

+ 

0.037C 

T 

1 

= 

-0. 8500 

31 

-0.0316 

♦ 

0,041 7 

I 

1 


-0.9500 

32 

-0,0296 


0.C465 

\ 

7 

= 

-1, 0590 
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33 

-C.0272 


0.0513 

I 

2 

- 

-1. 1500 

34 

-0.0243 

4 - 

0.0563 

I 

7 


-1.2500 

3b 

-C.0210 

4 

0.C61 3 

I 

7 

- 

-1.3500 

36 

-0.0171 

+ 

0.0664 

I 

1 

- 

-1.450 0 

37 

-0.0127 

+ 

0.0 715 

I 

1 


-1.5500 

38 

-C.0088 

4 

0.0966 

I 

1 

= 

-1.6625 

39 

0.0001 

4 

0.1046 

1 

J 

= 

-1, 7875 

40 

0.0103 

4 

0.112 5 

I 

1 


-1.9125 

41 

0.0218 

4 

0.1203 

I 

7 

= 

-2.0375 

42 

0.0346 

4 

0.1279 

I 

1 

= 

-2. 1625 

43 

0.0490 

4 

0.1352 

I 

1 

= 

-?,2H75 

44 

0.0649 

4 

0. 142 3 

I 

I 

= 

-2,4125 

45 

0.0825 

4 

0.1491 

I 

1 

= 

-2.5375 

46 

C. 1018 

4 

0,1555 

I 

f 

= 

-2.6625 

47 

0. 1232 

4 

0.1615 

I 

1 

= 

-2. 7875 

4d 

0. 1466 

4 

0.1671 

I 
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